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Abstract 

This  thesis  provides  a  thorough  review  of  conventional  matched-filter  radar 
theory  as  it  applies  to  aerosol  sensing  lidar.  Basic  matched-filter  radar  theory  and 
the  complex,  dense  aerosol  target  model  eventually  lead  to  a  general  derivation  of  the 
matched-filter  radar  receiver  response  to  a  dense  aerosol  target  environment.  The 
matched-filter  response  is  obtained  in  terms  of  a  two  dimensional  convolution  of  the 
target  scattering  function  and  the  radar  ambiguity  function.  The  range  and  radial 
velocity  resolution  of  various  radar  signals  were  compared  using  the  matched-filter 
radar  receiver  and  scattering  function  models.  Pulse  compression  radar  signals  were 
compared  to  the  simple  radar  pulses  currently  in  use  with  existing  lidar  systems, 
and  in  each  case,  the  pulse  compression  radar  signal  did  not  provide  a  significant 
improvement  in  combined  range  and  radial  velocity  resolution. 
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Pulse  Compression  for  Aerosol  Ranging  With  Coherent 
Pulse-Doppler  Lidar  Systems 


1.  Introduction 


1.1  Background 

Rapid  advances  in  coherent  laser  pulse-Doppler  radar,  commonly  known  as 
coherent  lidar,  have  resulted  in  numerous  wind  sensing  applications  such  as  mea¬ 
surement  of  wind  velocity  profiles  of  severe  storms,  clear  air  turbulence  (CAT),  wind 
shear,  and  aircraft  wake  vortices.  Some  of  these  applications,  such  as  aircraft  wake 
vortex  detection  and  wind  shear  detection,  require  range  resolutions  which  exceed 
the  capabilities  of  current  lidar  systems.  Extending  classical  pulse  compression  tech¬ 
niques  originally  developed  for  conventional  radio  frequency  radar  to  lidar  systems 
was  believed  to  offer  the  possibility  of  improving  the  combined  range  and  velocity 
resolution  in  comparison  to  current  lidar  systems  which  employ  simple  unmodulated 
radar  pulses. 

1.2  Problem 

Present  “state  of  the  art”  lidar  systems  for  wind  sensing  with  a  required  velocity 
resolution  of  approximately  1  m/s  are  limited  in  range  resolution  to  approximately 
200  m.  The  objective  of  this  thesis  was  to  research  classical  pulse  compression 
techniques  to  determine  their  potential  for  improving  the  combined  range  and  radial 
velocity  resolution  of  lidar  systems.  Four  basic  claisses  of  radar  signal  were  identified. 
The  resolution  in  combined  range  and  radial  velocity  for  each  class  of  radar  signal 
was  determined  in  terms  of  the  conventional  matched-filter  radar  receiver  response 
to  a  dense  aerosol  target  environment.  The  resolution  performance  for  each  class  of 
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pulse  compression  radar  signals  was  thoroughly  reviewed  and  compared  to  the  simple 
unmodulated  radar  pulse.  The  thesis  research  showed  that  for  a  dense  aerosol  target 
environment,  pulse  compression  radar  signals  do  not  provide  improved  resolution  in 
combined  range  and  radial  velocity  in  comparison  to  the  simple  unmodulated  radar 
signals  currently  in  use  with  existing  lidar  systems. 

1.3  Summary  of  Current  Knowledge 

Lidar  systems  have  become  a  powerful  tool  for  the  measurement  of  radial  wind 
velocities  with  numerous  wind  sensing  applications.  Extensive  research  efforts  which 
began  in  the  late  1960s  and  have  continued  to  the  present  have  resulted  in  lidar  sys¬ 
tems  with  CO2  continuous  wave  (CW)  transmitters  which  are  amplitude  modulated 
to  obtain  pulse  energies  up  to  IJ  with  pulse  repetition  frequencies  (PRF)  up  to  50 
Hz.  A  typical  system  operating  with  a  CO2  wavelength  of  10  /zm  and  with  a  pulse 
duration  of  1  fis  can  attain  a  radial  velocity  resolution  of  approximately  1  m/s  and 
a  radial  range  resolution  of  approximately  150-300  m  (24).  These  systems  employ 
optical  heterodyne  detection  which  requires  a  monochromatic  CW  leiser  source  which 
is  spatially  and  temporally  coherent.  The  coherence  time  of  the  Iciser  transmitter 
must  be  greater  than  the  coherence  time  of  the  received  signal  which  is  determined 
by  atmospheric  conditions.  The  atmosphere  is  modeled  as  a  diffuse  target  which 
consists  of  a  large  number  of  very  small  aerosol  particles  such  as  dust  and  water 
droplets.  When  a  pulse  of  coherent  laser  energy  is  transmitted  to  the  atmosphere, 
the  pulse  irradiates  a  volume  of  air  (defined  by  the  beam  width  and  pulse  duration) 
which  consists  of  a  large  number  of  these  aerosols.  According  to  Mie  theory,  the 
backscattered  radiation  of  the  received  signal  is  due  primarily  to  suspended  aerosols 
with  diameters  within  one  order  of  magnitude  of  the  laser  transmitter  wavelength 
(18).  For  example,  with  a  10  /xm  wavelength  transmitter,  the  backscattered,  radia¬ 
tion  is  due  to  aerosols  with  diameters  of  approximately  1-3  fim  (24,  18).  Doppler 
shift  measurements  of  the  backscattered  optical  signal  as  a  function  of  time  result  in 
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range-resolved  measurements  of  radial  (line  of  sight)  velocity. 

The  optical  signal  is  heterodyned  to  an  electrical  signal  where  conventional 
Doppler  signal  processors  originally  developed  for  millimecer-wave  radar  are  em¬ 
ployed.  The  received  optical  signal  backscattered  from  the  range  resolved  volume  of 
aerosol  particles  may  be  processed  to  obtain  the  velocity  mean  and  velocity  standard 
deviation  using  either  time  domain,  pulse-pair  processing  techniques  or  frequency  do¬ 
main,  fast  fourier  transform  (FFT)  techniques.  The  most  common  signal  processing 
approach  requires  a  complex  receiver  structure  of  in-phcise  and  quadrature  compo¬ 
nents  to  determine  the  correct  target  velocity  direction  relative  to  the  lidar  followed 
by  tr«nsient  digitizers  and  a  pulse-pair  or  spectral  processor. 

1.4  Assumptions 

A  considerable  degree  of  research  concerning  all  aispects  of  lidar  wind  sensing 
has  occurred  over  the  past  two  decades,  and  models  for  atmospheric  effects  and 
photodetector  quantum  noise  effects  have  been  rigorously  developed.  Atmospheric 
models  to  account  for  attenuation,  speckle,  and  refractive  turbulence  are  zissumed 
adequate  for  this  research  effort.  Considerable  agreement  found  in  the  literature 
indicates  the  validity  of  these  models.  The  net  effect  of  the  atmospheric  models  is 
a  received  optical  signal  which  is  statistically  modeled  as  a  time  varying  Rayleigh 
phasor  with  a  Rayleigh-distributed  amplitude  and  a  uniformly-distributed  phase. 
The  noise  current  of  the  quantum  noise  limited  photodetector  is  statistically  modeled 
as  additive  white  Gaussian  noise  over  the  narrow  bandwidth  of  the  detector,  and 
the  photodetector  output  current  due  to  the  received  optical  signal  is  modeled  by 
a  narrowband  Gaussian  random  process  (24).  Errors  which  may  arise  from  these 
statistical  models  or  simplifications  of  the  models  are  accounted  for  in  the  analysis 
as  necessary. 

In  addition  to  assumptions  related  to  the  atmospheric  and  photodetector  sta¬ 
tistical  models,  two  additional  assumptions  are  made.  The  first  assumption  is  based 
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on  Biernson  and  Lucy’s  comparison  between  coherent  and  direct  detection  lidar  per¬ 
formance  which  shows  that  for  low  signal  to  noise  ratio  {SNR  <  1)  coherent  (hetero¬ 
dyne)  detection  provides  superior  detector  output  SNR  performance  (4).  Therefore, 
since  low  SNR  is  prevalent  in  wind  sensing  applications,  the  analysis  will  be  based 
only  on  a  heterodyne  detection  model.  The  second  assumption  is  based  on  a  review 
of  current  lidar  transmitters  which  indicates  that  CO2  lidar  transmitters  operating 
at  9-12  fim  wavelengths  represent  the  most  mature  transmitter  technology  available 
for  lidar  today.  Therefore,  the  analysis  will  be  based  only  on  a  CO2  transmitter 
model  which  operates  with  a  wavelength  of  9-12  fim. 

1.5  General  Approach  and  Presentation 

Results  of  this  thesis  research  provide  a  thorough  review  of  conventional  matched- 
filter  radar  theory  as  it  applies  to  an  aerosol  sensing  lidar.  Chapter  two  provides  an 
historical  perspective  of  lidar  development  along  with  a  brief  discussion  of  current 
lidar  systems  and  signal  processing  techniques.  Chapter  3  provides  a  discussion  of 
the  basic  theory  of  conventional  matched-filter  radar  receivers,  the  complex  aerosol 
target  model,  and  the  matched-filter  radar  receiver  response  to  the  aerosol  target 
model.  In  Chapter  4,  each  of  the  four  basic  classes  of  radar  signal  are  evaluated  in 
terms  of  their  resolution  performance  in  a  dense  aerosol  target  environment. 
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11.  A  Review  of  Lidar  Developments 

Coherent  Lidar  systems  have  become  a  powerful  tool  for  the  measurement  of 
radial  wind  velocities  with  numerous  wind  sensing  applications.  This  chapter  pro¬ 
vides  a  basic  overview  of  coherent  lidar  technological  and  theoretical  developments 
occurring  over  the  past  two  decades.  The  chapter  begins  with  a  discussion  of  lidar 
historical  developments,  followed  by  a  simplified  lidar  system  description,  a  brief  dis¬ 
cussion  of  lidar  theory,  and  finally,  a  discussion  of  lidar  signal  processing  techniques. 

2.1  Historical  Development  of  Lidar  Systems 

In  1980,  Bilbro  of  the  NASA  Marshall  Space  Flight  Center  (MSFC)  provided 
a  detailed  overview  of  the  historical  developments  of  atmospheric  laser  Doppler  ve- 
locimetry  which  occurred  from  1964  to  1980  (6).  According  to  Bilbro,  laser  Doppler 
velocimetry  theory  baised  on  the  meeisurement  of  the  mean  frequency  of  a  backscat- 
tered  target  frequency  spectrum  Wcis  first  described  in  1964.  During  the  early  1970s, 
the  first  demonstrated  measurements  of  atmospheric  wind  velocity  occurred  using 
continuous  wave  (CW)  carbon  dioxide  CO2  laser  systems  with  operating  ranges  of 
approximately  200  m  (6,  22).  These  initial  CO2  CW  laser  systems  were  used  in  a 
variety  of  applications  in  wind  sensing  such  as  measuring  the  wind  velocity  profiles 
of  wind  shear,  dust  devils,  tornadoes,  and  aircraft  wake  vortices  (6,  22,  32,  17). 

During  the  1970s  and  soon  after  CW  systems  were  designed,  the  desire  to 
measure  clear  air  turbulence  (CAT)  motivated  the  development  of  coherent  pulsed- 
Doppler  lidar  systems  (6,  24).  These  cysterns  offered  a  tenfold  increase  in  range 
performance.  The  first  air-borne  pulsed  lidar  system  was  built  for  NASA/MSFC 
and  employed  a  master  oscillator  power  amplifier  (MOPA)  transmitter  configura¬ 
tion  capable  of  attaining  pulse  energies  up  to  30  mJ  and  operating  ranges  up  to  6 
km  (5,  7).  The  system  was  later  upgraded  with  the  capability  of  generating  two- 
dimensional  wind  fields  horizontal  to  the  aircraft.  The  technique  for  generating  these 
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two-dimensional  wind  profiles  was  based  on  the  measurements  of  wind  velocity  at  the 
same  location  from  two  directions.  This  upgraded  airborne  lidar  system  was  flown 
on  the  NASA  CV-990  during  several  studies  in  the  early  1980s  which  successfully 
measured  two-dimensional  wind  profiles  surrounding  severe  storms  (5,  7). 

According  to  Menzies  and  Hardesty,  the  development  of  pulsed  transversely- 
excited  atmospheric  pressure  (TEA)  CO2  Ijiser  transmitters  was  a  major  milestone 
of  Doppler  lidar  technological  development,  and  in  1981  the  National  Oceanic  At¬ 
mospheric  Administration  (NO  A  A)  began  using  the  first  pulsed  TEA  CO2  system 
to  measure  atmospheric  wind  velocities  (24).  This  TEA  laser  system  produced  a 
pulse  energy  of  100  mJ  with  a  pulse  repitition  frequency  (PRF)  of  10  Hz.  In  a 
continuing  effort  to  meet  the  pulse  energy  requirements  of  a  space-based  lidar  sys¬ 
tem,  the  NOAA  Wave  Propogation  Lab  (WPL)  later  upgraded  their  hybrid  TEA 
laser  transmitter  lidar  system  with  an  injection-locked  TEA  laser  capable  of  1  J 
pulse  energies  with  a  maximum  PRF  of  50  Hz  (13).  The  NOAA  system  with  the 
injection  controlled  TEA  transmitter  and  velocity  azimuth  display  (VAD)  scanning 
techniques  to  generate  three-dimensional  wind  velocity  profiles  with  color  graphics 
displays  represents  a  “state  of  the  art”  lidar  system  which  has  been  used  in  numerous 
wind  sensing  studies  (24). 

A  major  application  and  motivation  for  the  development  of  lidar  technology 
hais  been  the  concept  of  an  earth  orbiting  lidar  system  to  generate  globail  three- 
dimensional  wind  velocity  profiles  (24, 10).  Numerous  feasability  studies,  preliminary 
design  specifications,  and  advances  in  coherent  pulse  Doppler  lidar  technology  during 
the  past  decade  have  resulted  in  a  baseline  lidar  system  called  the  Laser  Atmospheric 
Wind  Sounder  (LAWS)  which  is  scheduled  to  become  part  of  NASA’s  Earth  Orbiting 
System  (EOS)  in  the  late  1990s  (23,  27,  40).  The  baiseline  LAWS  system  design 
specifications  include  a  CO2  laser  transmitter  operating  at  a  wavelength  of  9.11  fim 
with  3  (IS  pulse  duration,  10  J  pulse  energy,  and  a  PRF  of  up  to  10  Hz  (24,  10). 
More  recent  developments  in  solid-state  laser  transmitters  for  pulsed  lidar  may  result 
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in  a  solid-state  coherent  lidar  implementation  for  LAWS  (10,  2,  9).  For  example, 
Coherent  Technologies  Inc.  has  recently  developed  and  demonstrated  a  partially 
solid-state  Nd:YAG  laser  transmitter  design  (19).  Rapid  advances  in  lidar  technology 
have  resulted  in  accurate  and  reliable  coherent  pulsed-Doppler  lidar  systems  for 
three-dimensional  wind  field  measurements  with  the  eventual  goal  of  a  space-based 
lidar  fystem  generating  global  three-dimensional  wind  velocity  profiles. 

2.S  Lidar  System  Description 

Several  authors  have  described  the  theory  of  coherent  laser  pulsed-Doppler 
radar  developed  more  than  two  decades  ago  (4,  1).  The  simplified  block  diagram  of 
Figure  1  illustrates  a  typical  coherent  lidar  system  employing  a  MOPA  configuration. 
A  signal  trace  through  the  various  components  of  Figure  1  will  simplify  the  system 
description. 

Within  the  MOPA  configuration,  the  spatially  coherent,  monochromatic  CW 
optical  or  infrared  output  beam  is  amplitude  modulated  and  power  amplified  to  form 
an  optical  pulse  train  signal  with  a  carrier  frequency  equal  to  the  frequency  of  the 
master  oscillator.  The  optical  signal  is  transmitted  through  the  transmit/receive 
optical  duplexer  to  the  telescope  and  possibly  through  a  scanning  mirror  to  the 
target  which  in  this  application  is  the  aerosols  suspended  within  the  atmosphere. 
The  switching  electronics  of  the  optical  duplexer  establishes  the  time  measurements 
for  range  gating  and  protects  the  detector  from  the  high  transmission  intensities. 
In  a  monostatic  radar  configuration  as  illustrated  in  Figure  1,  the  same  telescope 
receives  the  optical  signal  backscattered  from  the  target.  The  velocity  of  the  target 
will  cause  a  shift  of  the  mean  frequency  of  the  received  signal  frequency  spectrum 
according  to  the  well  known  Doppler  equation: 


Fi  =  2Vrl\ 


(1) 
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Figure  1.  Simplified  Block  Diagram  of  a  Coherent  Pulsed-Doppler  Lidar 
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where  Fd  is  the  Doppler  frequency  shift,  K-  is  the  radial  velocity  of  a  target,  and  A  is 
the  wavelength  of  the  signal  transmitted.  Part  of  the  CW  master  oscillator  beam  is 
split-off,  frequency  translated,  and  optically  summed  with  the  received  target  echo 
signal,  and  the  optically  mixed  signal  is  focused  onto  the  photodetector  (4). 

The  photodetector  performs  as  a  square-law  device  with  negligible  second  order 
terms  converting  optical  energy  to  electrical  energy  and  creating  sum  and  difference 
frequency  terms  (4).  The  frequency  bandwidth  of  the  detector  limits  the  output  to 
only  the  frequency  difference  terms;  therefore,  the  output  of  the  photodetector  is  an 
electrical  signal  with  a  radio  frequency  (RF)  spectrum  and  a  mean  frequency  equal 
to  the  difference  between  the  received  signal  frequency  and  the  frequency  translated 
master  oscillator  frequency.  This  detection  procedure  is  commonly  known  as  optical 
heterodyne  detection. 

A  well  designed  coherent  optical  receiver  is  a  shot  noise  or  quantum  noise 
limited  receiver  in  which  the  primary  source  of  noise  is  the  shot  noise  generated 
by  the  mcister  oscillator.  The  additive  shot  noise  of  the  photodetector  is  modeled 
statistically  as  additive  white  Gaussian  noise  over  the  narrow  bandwidth  of  the 
detector.  With  a  monochromatic  CW  oscillator  of  sufficient  power  so  that  other 
noise  sources  such  as  dark  current  noise  and  optical  background  radiation  become 
negligible,  the  optical  receiver  can  achieve  an  average  signal  to  noise  ratio  expressed 
as: 


SNR  = 


yPr 

hvB 


(2) 


where 

SNR  =  average  signal  to  noise  ratio 
T}  =  photodetector  quantum  efficiency 
Pr  =  average  power  of  the  received  signal  with  the  same 
polarization  as  the  master  oscillator  signal  (W) 
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h  =  Planck’s  constant  (J-s) 

1/  =  optical  frequency  incident  on  photodetector  (Hz) 

B  =  photodetector  bandwidth  (Hz) 


According  to  Bachman  (1:11-16)  when  the  target  echo  signal  is  a  coherent, 
normally  incident  plane  wave,  the  output  of  the  photodetector  is  an  ensemble  average 
electrical  current  expressed  as: 

i{t)  =  ^{Pr^o  +  Pr{t)  +  2[PmoPrYf^COs[2x{Uo  "  fr)t  +  ^r(<)]}  (3) 


where 


1}  =  photodetector  quantum  efficiency 

q  =  electronic  charge  (C) 

Pmo  —  power  of  the  mcister  oscillator  CW  signal  (W) 

P,(<)  =  power  of  the  received  signal  with  the  same 

polarization  as  the  master  oscillator  signal  (W) 

=  frequency  of  master  oscillator  signal  plus  offset  (Hz) 
fr  =  frequency  of  received  signal  (Hz) 

=  phase  difference  between  the  received  signal 
and  the  master  oscillator  signal  (rads) 


The  second  term  of  Equation  (3)  is  the  frequency  modulated  signal  of  interest. 

After  heterodyning,  the  electrical  signal  is  input  to  the  signal  processor  and 
computer  where  either  time  domain  or  frequency  domain  techniques  originally  de¬ 
veloped  for  millimeter-wave  radar  are  employed.  Using  classical  time  domain  and 
frequency  domain  signal  processing  techniques,  target  information  such  as  range  and 
range-rate  is  obtained.  After  accounting  for  atmospheric  effects  and  the  receiver  SNR 
from  Equation  (2)  which  are  different  from  millimeter-wave  radar  receivers,  Bach- 
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man  concludes  that  the  statistical  problem  of  accurately  measuring  target  range  and 
target  radial  velocity  with  a  coherent  lidar  is  essentially  the  same  as  the  detection 
problem  using  conventional  millimeter-wave  radar. 


2.S  Lidar  Theory 

Conventional  millimeter-wavelength  and  laiser- wavelength  pulse-Doppler  radars 
measure  the  phase  difference  between  subsequently  transmitted  and  reflected  pulses 
to  determine  hard-target  range  and  radial  velocity;  whereas,  with  a  diffuse  target 
such  as  the  atmosphere,  Doppler  lidar  me2isures  the  Doppler  shift  of  each  received 
pulsed  defined  by  the  range-gate  of  the  system.  A  hard-target  is  simply  defined  as  a 
target  which  is  smooth  on  the  scale  of  the  laser  wavelength  transmitted  and  reflects 
the  transmitted  signal  without  amplitude  or  frequency  distortion.  A  diffuse  target 
is  rough  on  the  scale  of  the  transmitted  laser  wavelength  which  results  in  multiple 
reflections  of  the  transmitted  signal.  Several  authors  have  provided  the  radar  range 
equation  which  describes  the  hard-target  detection  scenario  (37, 1,  20)  For  the  hard- 
target  case  where  the  target  is  larger  than  the  radar  beamwidth  and  normal  to  the 
plane  of  the  optical  beam,  the  radar  range  equation  can  be  expressed  £is: 


Pr.{t)  =  Pt{t-2Rlc)p^r)seexp  -2  a{P!)dP! 


(4) 


where 

«■(<) 

m 

R 

c 

p 


=  expected  value  of  the  received  power  with  the 
same  polarization  as  the  master  oscillator  (W) 

=  range-delayed  power  profile  of  the  transmitted  signal  (W) 
=  range  to  the  target  (m) 

=  speed  of  light  (m/s) 

=  target  reflected  power  per  steradian  (sr)  divided  by  the 
incident  power  (l/sr) 
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Af  =  effective  area  of  the  antennae  aperture  (m^) 

J?je  =  optical  system  efficiency 

ot{R)  —  two-way  extinction  coefficient  (1/m) 

In  Equation  (4)  above,  some  simplifications  have  been  made.  The  received 
signal  power  Pr(<)  represents  an  average  power,  and  its  actual  value  is  slightly  less 
when  the  atmospheric  effects  of  speckle  pattern  fluctuation  and  refractive  turbulence 
are  included.  For  hard-target  detection  with  sufficient  SNR,  turbulence  effects  and 
the  aerosol  volume  backscatter  coefficient  ^(R)  (1/mxl/sr)  become  negligible.  The 
power  profile  of  the  return  signal  range  delayed  in  time  is  the  same  as  the  transn  itted 
power  profile;  however,  for  the  case  of  a  diffuse  target  such  as  a  volume  of  aerosols, 
the  returned  signal  power  profile  changes  significantly. 

For  the  diffuse  target  case,  the  radar  range  equation  (20:2)  can  be  expressed 
as: 

Pr{t)  =  Pt{t  -  2i?/c)/?(P)4|»/,r  exp  [-2  r  a(R')d^i  (5) 

Jc(t-T)/2  it  L  J 

where  t  is  the  pulse  duration.  Again,  the  effect  of  turbulence  has  been  ignored.  If 
strong  turbulence  is  present,  an  additional  turbulence  factor  to  reduce  Pr  is  required. 
The  ensemble  average  power  profile  Pr{t)  depends  on  the  cumulative  effect  of  the 
backscattered  radiation  from  a  large  number  of  aerosols  within  the  illuminated  vol¬ 
ume  which  is  defined  by  the  beam  width  and  the  pulse  duration.  The  integration 
over  R  of  Equation  (5)  indicates  that  at  time  t,  the  aerosols  within  a  slab  of  atmo¬ 
sphere  of  thickness  cr/2  and  centered  at  R  =  c{t/2  —  t/4)  backscatter  radiation  to 
form  the  power  profile  Pr{t)  (20).  Therefore,  the  range  resolution  defined  by  ct/2 
is  relatively  coarse  compared  to  a  conventional  pulse  detection  Doppler  radar.  For 
a  typical  CO2  transmitter  with  a  pulse  duration  of  1  jtzs,  the  range  resolution  is 
approximately  150  m.  If  the  pulse  duration  is  reduced  to  improve  range  resolution. 
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the  resultant  spread  of  the  frequency  spectrum  U  the  transmitted  signal  will  result 
in  a  significant  loss  of  frequency  (velocity)  resolution. 

Therefore,  the  classical  trade-off  between  range  and  velocity  resolution  as  de¬ 
termined  by  the  pulse  duration  of  the  transmitted  signal  exists  for  lidar  diffuse  target 
applications  2is  well  as  conventional  Doppler  radar  hard-target  applications.  Accord¬ 
ing  to  Menzies  and  Hardesty,  the  accuracy  of  Doppler  lidar  velocity  measurements 
depend  on  two  fundamental  parameters:  the  SNR  and  the  spectral  width  of  the 
backscattered  signal  (24). 

2.4  Lidar  Signal  Processing 

Although  this  thesis  does  not  directly  apply  lidar  signal  processing  techniques, 
a  brief  introduction  to  current  lidar  wind  sensing  signal  processing  techniques  is 
given  to  provide  a  complete  overview  of  lidar  theory.  Current  lidar  wind  sensing  sys¬ 
tems  typically  employ  three  stages  of  signal  processing:  amplification  and  complex 
demodulation,  analog  to  digital  (A/D)  conversion,  and  either  frequency  domain  or 
time  domain  digital  signal  analysis.  After  the  received  signal  from  the  photodetector 
is  amplified,  filtered,  converted  to  quadrature  components,  and  digitized,  the  signal 
parameters  of  interest  are  estimated  using  one  of  the  two  signal  processing  techniques 
developed.  These  two  techniques  are  Fourier  spectral  analysis  and  c(  ariance  analy¬ 
sis.  Fourier  spectral  analysis  is  a  relatively  straight  forward  method  of  obtaining  an 
estimate  of  the  power  spectral  density  (PSD)  of  the  received  signal  which  is  referred 
to  as  a  periodogram.  Three  parameters  of  the  periodogram  of  primary  interest  are 
the  following:  the  mean  power,  the  mean  frequency,  and  the  frequency  width.  In 
comparison  to  Fourier  spectral  analysis,  the  covariance  technique  is  a  less  compu¬ 
tationally  intense  method  of  obtaining  estimates  of  these  three  parameters.  In  this 
section,  complex  demodulation,  A/D  conversion,  and  PSD  parameter  estimation  for 
current  lidar  systems  are  briefly  reviewed. 
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Complex  demodulation  structures  such  as  the  one  illustrated  in  Figure  2  typi¬ 
cally  consist  of  a  bandpass  filter  broad  enough  to  receive  the  Doppler-shifted  signal 
Sri  {t)  from  the  photodetector  output.  The  signal  s,,  (t)  is  filtered  and  amplified  to 
produce  the  intermediate  signal  Srj{i).  The  signal  Sr2(i)  is  divided  into  two  chan¬ 
nels  called  the  I  and  Q  channels  for  in-phase  and  quadrature  respectively.  Mixing 
Srj(t)  with  the  single  frequency  terms  s\n{2irfot)  and  cos(2x/o<)  and  low  pass  fil¬ 
tering  forms  the  in-phase  and  quadrature  components  of  Sr^{t).  The  low  pass  filter 
(LPF)  of  each  channel  preserves  the  frequency  difference  signals  and  eliminates  the 
frequency  sum  signals  associated  with  the  mixing  process.  The  mixing  and  filtering 
of  Sr^it)  to  obtain  the  I  and  Q  components  of  Srj(<)  provides  a  means  of  determining 
the  direction  of  the  radial  velocity  relative  to  the  lidar.  After  sampling  the  I  and  Q 
channels,  the  digital  output  signal  can  be  written  as 

Srj  (nZ)  =  i{hTs)  -}-  q{nZ)  (6) 

where  T,  is  the  sampling  period  and  n  =  0, 1,  ...,iV  —  1.  The  digitized  output  signal 
Sr^{nTs)  is  sent  to  a  computer  with  the  software  routines  necessary  for  either  spectral 
or  covariance  processing  to  estimate  the  signal  parameters  of  interest. 

First  consider  the  case  of  Fourier  spectral  processing.  Discrete  Fourier  trans¬ 
form  (DFT)  methods  are  employed  to  obtain  the  periodogram  for  each  range  bin  of 
Sr3(nra).  The  width  of  the  range  bin  for  each  periodogram  derived  is  determined  by 
the  number  of  digital  samples  N  of  the  DFT.  Coherent  Technologies  Incorporated 
(CTI)  has  developed  a  Fourier  spectral  processor  for  lidar  wind  sensing.  The  spec¬ 
tral  signal  processing  of  the  CTI  lidar  (42:41-42)  can  be  summarized  in  the  following 
steps: 

Step  1.  Obtain  the  DFT  (nominal  length  N  =  128  samples)  of  Sr3(nT,)  expressed 
cis: 

N-i 

Z  (kfres)  =  ^r^inTs)  exp  {-j2‘KnTskfxes)  (7) 

n=0 
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Figure  2.  Complex  Demodulator 


where  /res  =  l/NT^  and  k  =  0,l,  ...,N  —  1. 

Step  2.  For  each  range  bin,  determine  the  periodogram  expressed  as 

S{kfres)  =  ^\Z(kfreBf  (8) 

Step  3.  Store  the  periodogram  for  each  range  bin  in  memory,  transmit  additional 
pulses  (perhaps  as  many  as  1000)  and  repeat  the  process  storing  all  periodogram 
data  in  memory. 

Step  4.  After  transmitting,  processing,  and  storing  the  data  during  the  time  interval 
in  which  the  random  process  describing  the  backscattered  signal  is  assumed  station¬ 
ary  (refer  to  Chapter  3  for  details),  the  periodogram  data  can  be  edited  and  averaged 
to  eliminate  spurious  signals.  The  edited  and  averaged  periodogram  for  each  range 
b.n  designated  as  5(/;/ies)  ^  l^rgc  number  of  averaged  pulses  is  expected  to  yield 

a  symmetric  Gaussian  spectrum.  Therefore,  the  central  moments  of  S{kfxes)  yield 
the  desired  signal  parameters  of  .ean  power,  mean  Doppler  frequency,  and  Doppler 
variance. 
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Step  5.  Compute  the  first  three  central  moments  of  S(fc/res)>  The  zero  moment 
estimate  is  expressed  as 

N-l 


S{kfxes) 


(9) 


k=0 


where  P{R)  represents  the  estimated  average  power  backscattered  from  range  bin  R. 
The  range  bin  R  is  defined  as  R=  {c/2){N  x  IT,).  The  estimate  of  the  first  moment 
of  5(fc/res)  converted  to  radial  velocity  coordinates  is  expressed  as 


1 


k 

27  P[R)T, 


{k  -  km)  s  [modAr(A:)]] 


(10) 


where  A  is  the  wavelength  transmitted  and  km  =  fcmaxj  the  spectral  coefficient  of 
the  greatest  magnitude.  Vr{R)  is  the  estimated  average  radial  velocity  at  a  given 
range  bin.  The  estimate  of  the  second  central  moment  or  the  variance  of  S{kfies) 
for  each  range  bin  is  expressed  cis 


^2  +  ^ 

iPim 


2 

S  [modAr(fc)] 


(11) 


The  estimate  of  the  velocity  width  is  given  as 


avlR)  = 


(12) 


The  first  three  moment  estimates  of  S{k /res)  designated  as  P{R),  Vr{R),  and  crv{R) 
provide  the  information  of  interest.  The  first  three  moments  provide  the  mean  power 
estimate  which  indicates  the  aerosol  density,  the  mean  aerosol  radial  velocity  estimate 
which  indicates  the  mean  radial  wind  velocity,  and  the  spectral  width  estimate  which 
indicates  the  degree  of  relative  wind  velocity  due  to  atmospheric  turbulence. 

Next  consider  the  covariance  processing  technique.  This  technique  is  analogous 
to  the  pulse-pair  signal  processor  implemented  in  conventional  microwave,  Doppler 
meteorological  radar.  Using  the  single-lag  autocovariance  processor  (SLA)  (11:91- 
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119)  one  obtains  the  first  three  estimated  moments  of  the  received  signal  PSD. 
The  instantaneous  power  backscattered  can  be  expressed  in  terms  of  the  squared 
components  of  Srj{nTs)  due  to  pulse  m  as 

PAnT,)  =  P{nT,)  +  Q\nZ)  (13) 

Averaging  the  instantaneous  samples  of  Pmi^Ts)  determines  the  mean  backscattered 
signal  power  expressed  as 

I  I  A/-1  N-l 

where 

P{R)  =  the  estimate  of  average  power  at  range  bin  R  (W) 

R  =  the  range  bin  width  where  R  =  {c/2){N  x  Tj) 

N  =  number  of  digital  samples  within  a  range  bin 
M  =  number  of  pulses  averaged 

The  average  signal  power  estimate  less  the  white  Gaussian  noise  term  N  of  the 
photodetector  can  be  expressed  as 

S{R)  =  P{R)  -  N  (15) 

From  Doviak  and  Zrnic  (11:91-119),  the  SLA  processor  forms  an  estimate  of  the 
complex  autocorrelation  function  expressed  as 

i  E  «r.(nT.)s„  {{n  +  l)r,)  (16) 

n=l 
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For  symmetric  spectra,  the  unbiased  mean  radial  velocity  estimate  can  be  expressed 
in  terms  of  the  argument  of  the  complex  autocorrelation  function  estimate  R{Ts)  as 

{r(T.)}  (17) 


The  spectral  width  by  of  an  assumed  symmetric  Gaussian  spectrum  can  be  expressed 
in  terms  of  ^(Tj)  and  S{R)  of  Equations  (15)  and  (16)  as 


In 


(18) 


where  the  sgn  term  identifies  possible  negative  terms  due  to  low  signal-to-noise  ratio. 
For  more  details  of  the  derivations  to  obtain  Equations  (17)  and  (18),  refer  to  the 
text  of  Doviak  and  Zrnic  (11:91-11). 

There  are  advantages  and  disadvantages  cissociated  with  each  signal  processing 
technique.  For  more  detailed  discussions  of  the  variances  associated  with  the  spectral 
and  covariance  estimators  and  comparisons  of  each  in  terms  of  their  advantages  and 
disadvantages,  refer  to  the  references  provided  (11,  42,  36). 

When  implementing  either  processing  technique,  the  lidar  range  resolution  is 
a  function  of  the  range  bin  length  R  with  N  samples,  the  pulse  duration  T,  and 
the  bandwidth  B  of  the  receiver  electronics.  For  example,  the  CTI  lidar  (42:41) 
operates  with  the  following  nominal  aerosol  mezisurement  parameters:  T  =  l/is, 
r,  =  10ns,  and  N  =  128  samples  centered  at  each  range  bin.  With  the  CTI  system, 
computing  the  first  three  moments  of  the  periodogram  (as  in  Steps  1-5)  provides  a 
range  resolved  profile  of  the  wind  the  lidar  pulse  encounters  eis  it  propogates  through 
the  atmosphere.  The  range  bin  consisting  of  128  samples  spaced  1.5  m  apart  in  range 
forms  a  range  bin  with  of  R=  {cj2){N  x  T,)  =  192m.  The  samples  taken  at  each 
range  bin  must  occur  within  the  coherence  time  of  the  received  signal;  otherwise, 
random  modulation  of  the  received  signal  broadens  the  frequency  spectrum.  This 


18 


aspect  of  the  backscattered  signal  is  discussed  in  more  detail  in  Chapter  3. 

Sampling  in  range  to  form  the  range  bin  reduces  the  range  r  'solution  of  the 
system  because  the  time  duration  required  to  form  the  N  sample  range  bin  R  results 
in  samples  from  a  pulse  which  hcis  propogated  in  range.  In  other  words,  during  the 
time  required  to  obtain  the  128  samples  within  a  range  bin,  the  pulse  effectively 
smears  through  additional  aerosol  scatterers  to  form  a  larger  scattering  volume.  The 
scattering  volume  is  defined  in  terms  of  the  lidar  beamwidth  and  the  spatial  extent 
of  the  transmitted  pulse  expressed  eis  cr/2.  In  addition,  the  bandwith  of  the  lidar 
will  determine  the  degree  of  spatial  averaging  of  the  scattering  volume.  In  other 
words,  a  finite  bandwidth  receiver  averages  the  backscattered  signal.  Therefore,  due 
to  pulse  width  and  bandwidth  constraints  of  the  lidar  system,  range  resolution  is 
typically  greater  than  the  range  bin  width  R  defined  by  the  DFT  length  N.  Range 
resolution  is  discussed  in  greater  detail  in  Chapters  3  and  4. 

This  section  has  provided  a  basic  overview  of  existing  lidar  signal  processing 
techniques  for  the  dense  aerosol  target  environment.  It  can  be  shown  that  Fourier 
spectral  processing  is  comparable  to  approximate  matched-filter  radar  receivers  when 
the  spectral  coefficients  of  the  DFT  are  equally  weighted  and  the  DFT  length  is 
comparable  to  the  time  duration  of  the  transmitted  rectangular  pulse  (14:82-84).  In 
the  next  chapter,  the  theory  of  matched-filter  radar  as  it  applies  to  a  dense  aerosol 
target  environment  is  developed,  and  matched-filter  signal  processing  is  compared 
to  the  Fourier  spectral  processing  technique. 
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III.  Theory  and  Methods 


3.1  Introduction 

The  primary  purpose  of  this  thesis  effort  is  to  consider  alternative  conventional 
radar  signals  to  the  simple  unmodulated  pulse  currently  in  use  with  existing  lidars 
to  determine  whether  these  alternative  signals  offer  some  means  of  improving  the 
resolution  of  aerosol  measurements.  To  accompll.h  this  purpose,  two  general  models 
are  required:  a  model  of  the  general  lidar  receiver  and  a  model  to  describe  the  aerosol 
target.  Section  3.2  discusses  the  basic  theory  of  matched-filter  (MF)  radar  receiver 
design.  Section  3.3  describes  resolution  properties  of  the  receiver  in  terms  of  the  radar 
ambiguity  function.  Section  3.4  characterizes  the  complex  target  which  consists  of  a 
dense  array  of  point  targets  distributed  in  range  and  frequency.  And  finally,  Section 
3.5  combines  these  two  models  to  describe  the  general  matched-filter  square-law 
receiver  (MFSLR)  response  to  the  complex  target.  The  developments  presented  in 
Sections  3.2  through  3.5  provide  the  essential  tools  necessary  for  determining  range 
and  velocity  resolution  of  various  radar  signals  with  respect  to  the  aerosol,  dense 
target  model. 

3.2  Matched-Filter  Theory 

Development  of  the  matched-filter  receiver  for  the  aerosol  sensing  lidar  was 
chosen  because  it  was  seen  to  provide  the  simplest  and  most  commonly  employed 
method  of  evaluating  various  signals  to  determine  their  measurement  capability  with 
respect  to  a  myriad  of  possible  target  scenarios.  Throughout  the  history  of  radar 
development,  the  ambiguity  function  bcised  on  the  matched-filter  response  hcis  been  a 
fundamental  measurement  tool  to  describe  radar  performance  in  terms  of  resolution, 
accuracy,  and  ambiguity.  The  matched-filter  receiver  is  designed  to  provide  optimal 
signal  detection  and  parameter  estimation  of  an  idealized  point  target  of  unknown 
range  and  radial  velocity  in  white  Gaussian  noise  as  illustrated  in  numerous  texts 
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(38,  26).  According  to  Van  Trees  (39:276-279)  the  matched-filter  receiver  followed 
by  a  square-law  envelope  detector  will  provide  the  maximum-likelihood  estimates  of 
point  target  range  and  radial  velocity  when  the  receiver  noise  is  Gaussian.  Whether 
the  matched-filter  receiver  followed  by  a  square-law  envelope  detector  is  the  optimal 
receiver  in  a  dense  target  environment  such  as  aerosols  is  less  clear  as  Rihaczek  and 
Van  Trees  indicate  (39,  30).  However,  a  basic  assumption  of  this  thesis  effort  is  to 
determine  resolution  performance  in  terms  of  a  conventional  matched-filter  radar 
receiver  design. 

Consider  the  simplified  radar  model  where  the  target  is  cissumed  to  be  a  sin¬ 
gle  point  target.  If  the  target  is  in  motion,  the  received  signal  is  represented  as  a 
time-delayed,  Doppler-shifted  version  of  the  transmitted  signal  corrupted  by  white 
Gaussian  noise  from  the  shot-noise  limited  photodetector.  For  this  model,  it  is  as¬ 
sumed  that  the  point  target  radial  velocity  remains  constant  over  the  pulse  duration. 
In  other  words,  higher  order  target  acceleration  terms  are  negligible.  (30:56-65)  Thus 
the  transmitted  pulse  is  expressed  a3 

s(<)  =  m(t)exp(j27r/ct)  (19) 


where 

s{t)  =  transmitted  signal  expressed  in  complex  form  (V) 
m{t)  =  complex  modulation  function  (V) 
fc  =  carrier  frequency  (Hz) 

The  use  of  complex  signal  notation  is  assumed  throughout  the  thesis,  so  a 
signal  is  assumed  to  be  complex  unless  stated  otherwise.  The  magnitude  of  m(t) 
is  the  envelope  of  the  real  signal.  The  transmitted  signal  represented  in  complex 
baseband  form  is  simply  expressed  cis 

s{t)  =  m{t)  (20) 
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Ignoring  changes  of  amplitude,  the  signal  backscattered  from  the  point  target  can 
be  represented  as 


$r{t)  =  m{t  -  Td)  exp[i27r(/c  +  fd){t  -  rd)]  +  n{t)  (21) 

where 

Td  =  cR/2 
fd  =  -2Vr/X 

and  where 

Td  =  range-delay  time  due  to  the  point  target  range  R  (s) 

fd  =  Doppler  frequency-shift  due  to  the  point 

target  radial  velocity  (Hz) 
c  =  speed  of  light  (m/s) 

A  =  transmitted  signal  wavelength  (m) 
n{t)  =  white  Gaussian  noise  signal  (V) 

Conversion  of  the  received  signal  to  complex  baseband  form  yields 

Sr{t)  =  m{t  -  Td)  exp  (j27r/di  +  j(f>)  +  n(i)  (22) 

where  ^  =  27r  (fc  -1-  fd)  rd  radians.  The  task  of  the  receiver  is  the  accurate  estimation 
of  the  parameters  Td  and  fd  cissuming  neither  parameter  is  known  apriori.  This  is 
the  classical  radar  problem  and  references  discussing  its  many  applications  abound 
(38,  30,  39,  26).  Radar  texts  commonly  demonstrate  the  optimality  of  matched-filter 
detection  of  the  backscattered  signal  from  a  single  point  target  in  white  Gaussian 
noise  (38:369-376).  If  the  target  is  stationary,  a  single,  linear  time-invariant  matched- 
filter  (MF)  is  adequate  to  achieve  maximum  likelihood  estimation  of  target  range¬ 
time  delay  r^.  The  MF  impulse  response  h{t)  which  maximizes  the  peak  signal  to 
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mean  noise  power  ratio,  can  be  expressed  in  terms  of  the  transmitted  signal  as 

h{t)  =  s*{-t)  (23) 

The  matched-filter  frequency  transfer  function  H{f)  is  expressed  through  the  Fourier 
transform  of  Equation  (23)  as 

m) = sv)  (24) 

First  consider  the  stationary  point  target  case  wher«.  fd  =  0,  and  the  only 
target  parameter  to  estimate  is  the  range  time-delay  T(J.  For  this  case,  Equation  (22) 
simplifies  to 

Sr{t)  =  m  (<  -  Td)  exp  {j (f>)  -b  n{t)  (25) 

where  ^  =  27r  fcTd.  One  can  easily  show  that  the  MF  response  ^(t)  can  be  expressed 
as 

/oo 

5r(t)s*(t  -  T)dt  (26) 

•OO 

The  form  of  Equation  (26)  shows  that  the  matched  filter  response  is  obtained  through 
a  cross-correlation  of  the  received  signal  with  a  time  shifted  replica  of  the  transmitted 
signal.  Substituting  Equation  (25)  into  Equation  (26)  and  neglecting  the  noise  term 
n{t)  yields 

y  (t)  =  exp  (;■  ^)  /  s{t-Td)s*{t-  r)  dt  =  exp  ( j  (f)R{T-  Td)  (27 ) 

J— OO 

where  R{t  —  Td)  represents  the  time-shifted  autocorrelation  function  of  the  trans¬ 
mitted  signal.  The  MF  response  can  also  be  expressed  in  terms  of  the  complex 
modulation  envelopes  as 

y(T)  =  exp  {j  (f>)  f  m{i-Td)m*  {t-  t)  dt  =  exp  (j  ^)R(t-  Td)  (28) 
Where  in  this  case  R(t  —  Td)  represents  the  time-shifted  autocorrelation  function  of 
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the  complex  modulation  envelope  of  the  transmitted  signal.  Since  the  phase  (f>  is 
unkown,  envelope  detection  is  required.  If  a  square-law  envelope  detector  follows  the 
MF,  the  matched-filter  square-law  receiver  (MFSLR)  output  can  be  expressed  as 

|!/(r)f  =  \R(t  -  nf  (29) 

where  the  peak  of  \R{t  —  Td)f  occurs  at  r  =  tj.  This  type  of  estimation  is  referred 
to  as  noncoherent  estimation  since  the  target  phase  information  is  eliminated. 

Now  consider  the  non-stationary  point  target  case  in  which  the  point  target 
parameters  rj,  /j,  and  ^  are  unknown.  For  this  case,  the  received  signal,  neglecting 
the  Gaussian  noise  term  can  be  expressed  cis 

Sr{t)  =  s{t  -  Td)  exp(i27r/dt  -f  jcf))  (30) 

Since  both  Td  and  fd  are  unknown,  an  array  of  matched-filters  with  incremen¬ 
tally  spaced  center-frequencies  as  illustrated  in  Figure  3  is  required  to  extract  the 
maximum-likelihood  estimates  of  Td  and  fd.  If  the  filters  of  the  MF  array  of  Figure 
3  are  spaced  with  infinitesmally  small  center-frequency  increments,  the  array  MF 
response  can  be  expressed  as 

y(T>/)=/  Sr{t)s*  {t-T)exp{-j2Trft)dt  (31) 

J—oo 


The  form  of  Equation  (31)  shows  that  the  MF  response  consists  of  the  cross¬ 
correlation  of  the  received  signal  with  a  time-shifted,  frequency-shifted  version  of 
the  transmitted  signal  (30:69).  The  received  signal  of  Equation  (30)  substituted  into 
Equation  (31)  yields 

y{T,  f)  =  exp  {j(l>)  I  s{t-  Td)  s*  {t  -  t)  exp  [-;27r  (/  -  fd)]  dt  (32) 

— OO 
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Figure  3.  Matched-Filter,  Square-Law  Receiver 
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Since  envelope  detection  is  typically  employed,  the  phase  term  is  neglected,  and  the 
MF  response  can  be  expressed  as 

/oo 

5  (t)  s*  (<  -  (t  -  Td)]  exp  [-j27r  (/  -  fd)]  dt  (33) 

-OO 

With  the  change  of  variables 

t'  =  T  —  Td  (34) 

/  =  S-U  (35) 

the  MF  response  of  Equation  (33)  simplifies  to 

y{T\  f')  =  J  (0  (“i27r(/)  dt  (36) 

The  variable  translations  of  Equations  (34)  and  (35)  shift  the  MF  response  from  the 
target  location  {rdjfd)  in  the  (r,/)  plane  to  the  origin  of  the  (t\/)  plane. 

Equation  (36)  represents  the  MF  response  neglecting  the  unknown  phase  term 
and  the  noise  term  of  Equation  (25).  This  function,  first  introduced  to  radar  theory 
by  Woodward  (41:118-119),  is  called  the  signal  ambiguity  function.  For  a  complex 
transmitted  signal  s(t),  the  ambiguity  function  is  defined  eis 

/oo 

-  t)  exp{-j2irft)dt  (37) 

-OO 

where  the  prime  notation  has  been  dropped.  The  ambiguity  function  represents  the 
MF  response  with  its  axes  centered  at  the  target  location  {Tdyfd)-  In  terms  of  the 
complex  low  pass  envelope  of  s(t),  the  ambiguity  function  is  expressed  as 

/CO 

m{t)m*{t  -  t)  exp{-j2Trft)dt  (38) 

-OO 
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The  ambiguity  function  can  also  be  expressed  in  its  symmetrical  form  as 


/OO 

m{t  +  Tl2)m*{t  —  T  l2)exp{—j2Trft)dt  (39) 

-00 


The  MFSLR  response  (as  in  Figure  3)  is  the  squared-modulus  of  x('^>/) 
designated  with  the  symbol  0(t,  /)  where 


^(’■j  /)  =  /  ”^(^  +  'r/2)m*(i  —  Tf2)exp{—j2'!rft)dt 

*/— 00 


(40) 


0(t,/)  is  also  commonly  referred  to  as  the  ambiguity  funtion,  and  the  three  dimen¬ 
sional  plot  of  the  surface  of  ^(t,/)  above  the  (r,/)  plane  is  called  the  ambiguity 
diagram. 

Several  mathematical  properties  of  the  ambiguity  function  0(r,  /)  can  be  found 
in  numerous  radar  texts  (39,  38).  These  properties  are  typically  given  in  terms  of  a 
transmitted  signal  ${t)  which  is  normalized  with  respect  to  its  energy  such  that 


dt  =  1 


(41) 


If  the  transmitted  signal  is  normalized  such  that  Equation  (41)  is  satisfied,  some 
useful  properties  of  the  ambiguity  function  (39,  38)  can  be  stated  in  the  following 
equations: 


Ambiguity  Volume 
Max{0(r,/)} 
0(t,O) 

^(0,/) 


0(0,0)  =  1 

roo  2 

/  m(<)m*(<  —  T)di 

4/ -00 

/OO 

m^(t)  exp(— j27r  }t)di 

•OO 

f) 


OO 

0(— T 


(42) 

(43) 

(44) 

(45) 

(46) 


Equation  (42)  is  known  as  the  volume  invariance  property  which  indicates  that  the 
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volume  of  a  normalized  ambiguity  function  always  equals  unity.  Equation  (43)  indi¬ 
cates  that  the  maximum  of  the  normalized  ambiguity  function  always  occurs  at  the 
origin  and  equals  unity.  Equations  (44)  and  (45)  indicate  the  form  of  the  ambigu¬ 
ity  function  across  the  range-time  delay  and  frequency  axes  respectively.  Equation 
(46)  is  a  symmetry  relationship.  The  volume  invariance  property  and  the  maxi¬ 
mum  amplitude  property  describe  the  fundamental  limitation  of  radar  resolution 
performance. 

Another  useful  function  in  radar  theory  is  the  cross  ambiguity  function  (26:286) 
defined  as 

/OO 

Si{t)sl{t  -  t)  exp{-j27rft)dt  (47) 

-<X> 

Substituting  Sr{t)  as  expressed  in  Equation  (25)  for  si{t)  and  substituting  the  trans¬ 
mitted  signal  for  S2{t),  the  cross  ambiguity  function  represents  the  MF  response 
given  as 

Xc(r,  /)  =  exp{j(f>)x  (r,  /)  +  Xn(r,  /)  (48) 

where 

Xn{rJ)=  [  n{i)s"{t-T)exp{-j2nft)dt  (49) 

Equation  (48)  represents  the  MF  response  in  terms  of  the  signal  component 
and  Gaussian  noise  component  of  the  received  signal.  The  cross  ambiguity  function 
represents  a  statistic  for  estimating  and  fd-  If  fd  is  known,  the  peak  of  \Xc{T^fd)\^ 
along  the  t  axis  determines  t^.  If  td  is  known,  the  peak  of  \XciTd,f)f  along  the 
/  axis  determines  fd-  If  neither  rj  nor  fd  are  known,  the  peak  of  \Xc{T,f)f  in  the 
(r, /)  plane  determines  rj  and  fd  with  maximum  likelihood  (8:8-58-8-59). 

3.3  Radar  Resolution 

Radar  resolution  is  defined  in  terms  of  the  radar’s  capability  to  separate  signals 
due  to  point  targets  which  are  separated  in  range  or  in  radial  velocity  or  in  both  range 
and  radial  velocity.  Radar  resolution  is  analyzed  baised  on  the  premise  that  the  signal 
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can  be  detected  in  the  receiver  noise;  otherwise,  analysis  of  resolution  capability 
is  irrelevant.  In  general,  the  signal  ambiguity  function  without  considering  noise 
effects  is  required  for  determining  radar  resolution  performance  (30:1-9).  Recall  from 
Equations  (34)  and  (35)  of  Section  3.2  that  the  definition  of  the  ambiguity  function 
x(t, /)  is  in  terms  of  the  MF  response  translated  to  the  target  coordinates  (rj, /j). 
Substitution  of  Equations  (34)  and  (35)  into  Equation  (38)  yields  the  expression  for 
the  MF  response  in  the  (t,  /)  plane  which  can  be  expressed  as 

y(T,  /)  =  Qxp{j(f>)x{r  -TdJ-  fd)  (50) 

Next,  consider  the  case  where  numerous  point  targets  are  imaged.  In  a  dense 
target  environment,  the  radar  beam  illuminates  many  closely-spaced  point  targets. 
If  cr{Tdnifdn)  represents  the  backscattering  cross  section  of  the  nth  point  target  and 
<f>{Td„,fd„)  represents  the  phase  of  the  nth  point  target,  the  MF  response,  neglecting 
the  Gaussian  noise  term,  becomes  a  superposition  of  the  MF  responses  from  the 
backscattered  signals  from  each  point  target.  The  resultant  MF  output  can  be 
expressed  as 


N 


!/(t, /)  =  S  /j„) exp  /*)]  X  (t  -  Ti,, /  -  /j.) 


(51) 


n=l 


where 

N  =  number  of  particles  illuminated 

(T^^^{Tdnifdn)  =  reflectivity  of  the  nth  point  target 

<t>{Tdni  fdn)  —  unknown  phase  of  the  nth  point  target  (rads) 

In  the  limit  cis  the  number  of  particles  becomes  infinite,  the  target  becomes  contin¬ 
uous  in  range  and  radial  velocity,  and  the  MF  response  is  given  as 

y{r,  =  J  J  ^  Id)  exp  \j(i>{rd,  fd)]  x{t  -  nj  -  fd)dTddfd  v52) 
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The  form  of  Equations  (51)  and  (52)  shows  that  the  MF  response  is  expressed  in 
terms  of  a  coherent  superposition  of  the  amplitudes  and  phases  of  each  point  target. 
As  a  result,  the  MF  response  as  expressed  in  Equations  (51)  and  (52)  is  called 
coherent  imaging  (8:7-12-7-13).  As  Rihaczek  explains,  in  a  dense  target  environment 
where  the  radar  illuminates  numerous  point  targets,  the  differential  ranges  between 
point  targets  are  essentially  random.  Therefore  attempts  to  determine  the  exact 
coherent  superposition  of  amplitudes  and  phases  to  describe  the  MF  response  is 
unrealistic  (30:111-114). 

To  proceed  with  the  discussion  of  resolution,  a  statistical  approach  is  required. 
Rihaczek  shows  that  the  dense  target  environment  requires  a  statistical  description 
of  the  target  backscatter  and  its  corresponding  MF  response  in  terms  of  the  mean 
target  backscatter  cross  section  a  and  the  mean  power  response  of  the  MF  y^.  The 
mean  power  response  can  be  obtained  with  the  MFSLR  illustrated  in  Figure  3.  If 
o{Td„,fd^)  represents  the  average  cross  section  of  the  nth  target  where  Td„  and  fd„ 
represent  the  mean  range  time-delay  and  Doppler  frequency-shift  of  the  nth  target, 
Rihaczek  shows  that  the  MFSLR  response  can  be  expressed  as 

\yi'rJ)f  ^  ^  (“T  -  -  fd„)  (53) 

n=l 

Again,  in  the  limit  eis  the  point  targets  become  continuous  in  range  and  frequency, 
Rihaczek  expresses  the  average  MFSLR  response  cis 

/)l^  -  /  -uj-  IdWddJ^  (54) 

Equation  (54)  describes  a  two-dimensional  convolution  of  the  average  scattering 
function  with  the  waveform  ambiguity  function.  Unfortunately,  limitations  exist 
concerning  the  validity  of  Equations  (53)  and  (54)  in  terms  of  imaging  aerosols. 
In  the  Rihaczek  derivation  of  these  equations,  he  states  the  assumption  that  the 
target  consists  of  a  fixed  distribution  of  scatterers  in  range  and  radial  velocity.  In 
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other  words,  the  orientation  of  the  particles  with  respect  to  the  radar  does  not 
change  significantly  over  the  time  duration  of  the  transmitted  pulse.  (30:335-338) 
This  limitation  is  especially  bothersome  when  imaging  aerosols  which  rapidly  change 
their  orientation.  If  the  pulse  duration  of  the  transmitted  signal  is  short  enough  (i.e. 
within  the  coherence  time  «  1  —  10/xs)  such  that  significant  changes  of  the  aerosol 
orientations  with  respect  to  each  other  do  not  occur  during  the  meaisurement  interval, 
Equations  (53)  and  (54)  are  valid  approximations  of  the  mean  MFSLR  response. 
Later  in  Sections  3.4  and  3.5  a  similar  but  more  general  result  for  the  mean  power 
MFSLR  response  is  obtained.  For  now,  eissuming  Equations  (53)  and  (54)  are  valid, 
they  can  be  used  to  illustrate  some  useful  concepts  of  resolution. 

The  form  of  Equation  (54)  shows  that  the  average  MFSLR  response  is  obtained 
through  a  superposition  of  0(r,  /)  for  each  point  target.  Also,  Equation  (54)  pro¬ 
vides  a  simple  analytical  description  of  the  mean  MFSLR  response  in  a  dense  target 
environment.  One  can  think  of  the  average  scattering  function  ^(rj,  f^)  as  a  two- 
dimensional  input,  the  ambiguity  function  6{t,  f)  as  an  impulse  response  of  a  two- 
dimensional  linear  filter,  and  the  mean  MFSLR  response  as  the  two-dimensional 
output  of  the  filter.  If  the  ambiguity  function  could  be  represented  as  the  tv'o- 
dimensional  impulse  function  S{T,f)  its  substitution  into  Equation  (54)  would  yield 


\y{T,  f)f  =  (55) 

The  two-dimensional  delta  ambiguity  function  represents  a  radar  with  infinitesimal 
resolution  in  time  and  frequency  or  alternatively  in  range  and  radial  velocity.  Based 
on  this  idealized  response,  a  radar  meteorologist  can  infer  common  parameters  of 
interest  pertaining  to  the  target  such  as  backscatter  density,  mean  velocity,  and 
velocity  spread.  These  parameters  are  obtained  from  tlie  complete  knowledge  of 
a(T,/).  It  is  well  known  that  the  delta  ambiguity  function  is  unrealizable  in  both 
the  T  and  /  axes  simultaneously.  A  realizable  radar  ambiguity  function  will  have 
finite  widths  At  and  A/  in  range  time-delay  (RTD)  and  Doppler  frequency-shift 
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(DFS).  It  is  clear  that  target  resolution  is  only  as  good  cis  the  filter  resolution  with 
which  it  is  observed;  therefore,  the  shape  of  the  ambiguity  function  0(t,  /)  and  its 
distribution  over  the  (r,  /)  axis  is  fundamental  to  determining  resolution. 

Three  commonly  referenced  resolution  criteria  are  the  Woodward,  the  Rayleigh, 
and  the  half-power  criteria.  The  Woodward  resolution  criterion  (41:118-119)  for  RTD 
resolution  width  can  be  expressed  as 

0(0,0) 

Equation  (56)  represents  the  central  width  Ar  in  which  the  energy  of  the  ambiguity 
function  is  located  with  respect  to  the  t  axis.  Assuming  the  transmitted  waveform 
is  normalized  with  respect  to  its  energy  and  substituting  Equation  (42)  of  Section 
3.2  into  Equation  (56)  yields 


At  =  /  0(t,  0)dr 

J  —  OO 

Substituting  Equation  (44)  of  Section  3.2  into  Equation  (57)  yields 


(87) 


—  T)dt 


dr 


(88) 


If  M{f)  represents  the  Fourier  transformation  of  m(t).  At  of  Equation  (58)  can  be 
expressed  as 

A 

mf)\'<if  (59) 

-00 

As  an  example,  let  the  normalized  complex  envelope  of  the  transmitted  signal 
spectrum  be  expressed  as 

Mf)\  =  (60) 
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where  B  is  the  bandwidth  and  where  the  function  vect{f/B)  is  defined  as 


rect 


1  i/i<f 

0  elsewhere 


(61) 


Substituting  Equation  (60)  into  Equation  (59)  yields 


df  =  l/B 


(62) 


Therefore,  the  total  resolution  width  in  RTD  is  inversely  proportional  to  the  signal 
bandwidth  B.  As  the  signal  bandwidth  increjises,  the  RTD  resolution  improves 
proportionately.  Similarly,  the  Woodward  DFS  resolution  criterion  for  a  normalized 
ambiguity  function  is  defined  as 


A/=  r  e{oj)df 

CO 


(63) 


Substituting  Equation  (45)  from  Section  3,2  into  Equation  (63)  and  again  assuming 
the  signal  energy  is  normalized,  DFS  resolution  can  be  written  as 


A/=  /  dt  (64) 

«/— OO 

And  a  similar  example  for  DFS  resolution  width  follows  where  the  transmitted  signal 
envelope  is  expressed  as 

=  ^rect(t/r)  (85) 

where  m{t)  is  a  normalized  rectangular  pulse  of  time-duration  T.  The  DFS  resolution 
width  is  obtained  from  the  following  equation: 

A/  =  |5^rect(l/r)|‘  <((  =  1/r  (66) 
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full  width  across  the  r  or  /  axis  respectively  such  that  0{T,f)  =  1/2.  The  RTD 
resolution  width  is  found  by  solving  for  At  in  the  following  equation: 


(71) 


Similarly,  the  DFS  resolution  width  is  found  by  solving  for  A/  in  the  following 
equation: 


e 


(72) 


For  the  simple  rectangular  pulse,  the  RTD  resolution  width  and  the  DFS  res¬ 
olution  width  is  usually  approximated  respectively  as 


At  =  Tc^l/B  (73) 

A/  =  1/Tc^B  (74) 


In  general,  range  delay  resolution  width  is  inversely  proportional  to  the  signal  band¬ 
width  and  frequency  resolution  width  is  inversely  proportional  to  the  signal  time 
duration  (30:93-98).  One  can  define  a  nominal  resolution  cell  with  dimensions  in 
frequency  and  range-time  delay  as 


At  A/  =  (75) 

where  TB  is  referred  to  in  radar  literature  as  the  time-bandwidth  product  of  the 
transmitted  waveform.  The  RTD  resolution  expressed  in  range  coordinates  is 

AR=cx  At/2  (76) 

And  the  conversion  of  DFS  resolution  to  radial  velocity  resolution  becomes 

AK  =  A  X  A//2  (77) 
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Substitution  of  Equations  (76)  and  (77)  into  Equation  (75)  yields  the  nominal  reso¬ 
lution  cell  with  dimensions  in  range  and  radial  velocity  expressed  as 

AHAK  =  7X3^  (’8) 

The  form  of  Equation  (78)  indicates  that  the  area  of  the  nominal  resolution  cell 
is  directly  proportional  to  the  wavelength  transmitted  and  inversely  proportional  to 
the  signal  time-bandwidth  product.  The  resolution  cell  definition  expressed  in  either 
Equation  (75)  or  in  Equation  (78)  provides  a  useful  analysis  tool  when  determining 
waveform  resolution  performance;  however,  the  resolution  cell  definition  alone  does 
not  provide  a  complete  description  of  resolution  performance.  The  resolution  mea¬ 
sures  At  and  A/  are  based  on  independent  measurements  of  range-time  delay  and 
frequency  of  the  central  lobe  of  the  ambiguity  function.  In  other  words,  Ar  is  the 
resolution  width  in  t  for  a  known  DFS  fj,  and  A/  is  the  resolution  width  in  /  for  a 
known  RTD  t^.  As  Rihaczek  (30:93-99)  and  Woodward  (41:118-124)  explain,  when 
combining  the  two  resolution  measurements  to  form  the  nominal  resolution  cell  one 
may  conclude  that  any  waveform  with  large  time-bandwidth  product  TjB  >  1  will 
result  in  a  smaller  resolution-cell.  However,  due  to  the  amplitude  constraint  of  Equa¬ 
tion  (43)  of  the  ambiguity  function  and  the  volume  constraint  of  Equation  (42)  of 
the  ambiguity  function,  waveforms  of  large  time-bandwidth  products  will  necessarily 
exhibit  ambiguities  elsewhere  in  the  (r,  /)  plane  which  can  mask  the  target  response 
of  the  resolution  cell  centered  at  the  origin  of  the  (r,  /)  plane. 

Rihaczek  introduces  the  term  self-clutter  to  describe  the  interference  of  un¬ 
wanted  targets  due  to  waveforms  with  ambiguity  functions  which  have  significant 
side-lobes  or  additional  ambiguities  in  the  form  of  near  replicas  of  the  central  spike 
in  the  (r,/)  plane  (30:92)  Therefore,  radar  resolution  in  a  dense  target  environment 
must  be  determined  in  terms  of  the  resolution  cell  definition  and  the  additional  ambi¬ 
guities  described  as  self-clutter  which  can  possibly  mask  the  resolution-cell  response. 


36 


In  this  section,  it  has  been  shown  that  the  MF  response  in  a  dense  target  en¬ 
vironment  is  described  statistically  in  terms  of  the  target’s  mean  backscatter  cross 
section  a(rd,  /^)  and  the  corresponding  mean  power  MFSLR  response.  This  statis¬ 
tical  average  approximation  of  the  target  backscatter  and  the  MFLSR  response  is 
valid  with  the  condition  that  the  point  targets  do  not  change  their  orientation  sig¬ 
nificantly  over  the  pulse  duration.  Since  the  mean  power  MFSLR  response  as  given 
in  Equation  (54)  describes  a  two-dimensional  convolution  of  the  scattering  function 
CT(rd,/j)  with  the  ambiguity  function  resolution  in  the  average  RTD  and 

average  DFS  is  defined  in  terms  of  the  central  spike  widths  At  and  A/  of  0(t, /). 
Three  radar  resolution  criteria  which  measure  At  and  A/  of  0(t,  /)  were  given.  In 
general,  target  resolution  will  require  knowledge  of  both  the  signal  ambiguity  func¬ 
tion  and  the  statistical  characteristics  of  the  target.  The  next  section  provides  a 
more  detailed  discussion  of  the  statistical  characteristics  of  the  aerosol  target. 

3.4  Target  Model 

The  discsussion  of  Section  3.3  led  to  a  statistical  model  of  the  dense  target 
environment  in  terms  of  the  target’s  mean  cross  section  This  statistical 

model  and  the  corresponding  MFSLR  response  is  valid  as  long  as  the  point  target  dis¬ 
tribution  in  Td  and  fd  which  composes  the  dense  target  does  not  change  significantly 
over  the  pulse  duration.  In  other  words,  changes  in  the  point  scatterer  distribution 
in  Td  and  fd  must  be  small  enough  that  oijdi  Id)  a  meaningful  quantity  to  describe 
the  target.  Since  aerosol  particles  rapidly  change  their  orientation  on  the  scale  of 
the  transmitted  laser  wavelength,  this  model  is  valid  only  for  short  duration  pulses 
(i.e.  fts  Ifis).  To  consider  the  more  general  case  where  the  orientation  of  the  aerosols 
changes  significantly  over  the  pulse  duration,  a  more  rigorous  statistical  model  of 
the  dense  aerosol  target  is  required. 

As  in  Section  3.3,  the  aerosol  target  model  can  be  described  as  a  randomly  dis¬ 
tributed  array  of  point  targets  which  are  distributed  in  range  and  velocity.  Aerosol 
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particles  of  a  diameter  on  the  order  of  the  optical  signal  wavelength  transmitted 
will  backscatter  most  of  the  energy  received.  Although  some  degree  of  molecular 
backscattering  occurs,  it  is  generally  assumed  to  be  negligible  in  comparison  to  the 
aerosol  backscattering  (18:67).  For  example,  under  typical  clear-air  atmospheric 
conditions,  most  of  the  energy  received  from  a  lidar  transmitting  at  the  10/im  wave¬ 
length  will  be  backscattered  from  aerosols  of  diameter  1 — 3/im  (24:450).  The  number 
density  of  these  particles  within  the  spatial  extent  of  the  pulse  transmitted,  depends 
on  a  multitude  of  factors  such  as  geographic  location,  altitude,  weather  conditions, 
season,  etc.,  and  will  typically  include  several  thousand  aerosols  (25).  From  a  the¬ 
oretical  point  of  view,  one  can  generally  assume  that  the  lidar  beam  will  illuminate 
an  infinite  number  of  point  targets  with  a  nominal  separation  on  the  order  of  a 
wavelength  for  lO^m  lidar  (15). 

From  Goodman  (12:347-356)  one  finds  that  when  monochromatic  light  illumi¬ 
nates  a  rough  surface  (rough  on  the  scale  of  the  wavelength  A)  the  reflected  image 
has  a  granular  appearance  due  to  numerous  bright  and  dark  spots.  This  effect  is 
commonly  referred  to  as  laser  speckle.  The  randomly  placed  aerosol  particles  within 
the  volume  which  the  lidar  beam  illuminates  individually  radiate  incoherent  fields  of 
random  phases.  At  the  receiver,  the  superposition  of  the  fields  backscattered  from 
each  of  the  aerosols  within  an  illumination  volume  will  form  the  resultant  speckle 
pattern.  The  illumination  volume  of  the  aerosols  is  defined  by  the  beam  width  and 
the  spatial  extent  of  the  transmitted  pulse  equal  to  cr/2  where  T  represents  the 
time  duration  of  the  pulse  transmitted.  The  received  signal  from  a  dense  target  such 
as  aerosols  can  be  expressed  as 

CO 

Sr{t)  =  exp  {j2Trfdj)  (79) 

n=l 


where 

Sr{t)  =  complex  received  signal  at  baseband  (V) 
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<7^^  =  reflectivity  of  the  nth  point  target 

4>n  =  unknown  phase  of  the  nth  point  target  (rads) 

m{t  —  Tj)  =  complex  modulation  signal  transmitted  (V) 

=  range-time  delay  of  the  nth  point  target  (s) 
fci„  =  Doppler  frequency  shift  of  the  nth  point  target  (Hz) 

From  Hardesty  (14:58-60),  the  backscattered  signal  Sr(t)  as  expressed  in  Equation 
(79)  represents  a  superposition  of  the  backscattered  signals  from  each  aerosol  point 
target.  The  phase  of  each  aerosol  is  statistically  described  by  a  uniform  distribu¬ 
tion  with  the  probability  density  function 

^  if  -TT  <  <  TT 

0  otherwise 

Therefore,  the  backscattered  signal  from  each  aerosol  point  target  has  a  random,  uni¬ 
formly  distributed  phase  and  the  backscattered  signal  Sr(i)  expressed  in  Equation 
(79)  represents  a  superposition  of  the  backscattered  signals  of  random  pheise  from 
each  point  target.  As  a  result,  for  a  given  time  ti  after  signal  transmission,  the 
received  signal  Sr(<i)  is  described  statistically  in  terms  of  its  amplitude  and  phase. 
The  amplitude  of  the  resultant  signal  \sr{t)\  is  described  statistically  as  a  Rayleigh 
random  variable  (14:58-60)  with  the  probability  density  function 


(SO) 


P{(f>n)  =  < 


Pi\Sr{h)\) 


|5r(<l)l  >  0 

\Sr{tl)\  <  0 


(81) 


where  is  the  mean  squared  value  of  |sr(ii)|.  The  resultant  pheise  of  Sr(ti)  desig¬ 
nated  as  <f>r{ii)  is  described  statistically  as  a  uniformly  distributed  random  variable 
with  the  uniform  probability  density  expressed  in  Equation  (80). 

As  Hardesty  and  Menzies  explain,  if  a  lidar  of  a  single  carrier  frequency  fc 
illuminates  a  volume  of  dense,  randomly  distributed  aerosols  which  are  stationary 
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with  respect  to  each  other,  the  resultant  backseat tered  signal  Sr{t)  will  consist  of  a 
sinusoidal  signal  of  the  carrier  frequency  fc  plus  or  miniis  the  Doppler  frequency-shift 
with  a  constant  Rayleigh  distributed  amplitude  and  constant  uniformly  distributed 
phase.  Under  real  atmospheric  conditions,  the  turbulence  and  wind  shear  within  the 
scattering  volume  cause  the  aerosols  to  move  relative  to  each  other.  The  relative 
aerosol  motion  within  the  scattering  volume  causes  the  aerosol  particles  to  change 
their  orientation  with  respect  to  the  lidar  receiver.  As  a  result,  the  superposition 
of  the  backscattered  signals  change  with  time  causing  the  random  amplitude  and 
random  phase  of  Sr{t)  to  also  change  with  time.  Therefore,  Sr(<)  is  described  as 
a  random  process.  The  changes  in  random  amplitude  and  random  phase  of  Sr{t) 
effectively  broaden  the  bandwidth  of  Sr{t)  (24:453). 

The  temporal  correlation  function  of  Sryt)  expressed  as 

Rsr{r)  =  E  [sr{t)s;{t  +  t)]  (82) 

provides  a  measure  of  how  rapidly  the  received  signal  changes  with  time.  If  the 
received  signal  is  assumed  to  represent  a  stationary  random  process  (this  representa¬ 
tion  of  Sr{t)  is  discussed  in  more  detail  in  the  following  paragraphs)  then  the  Fourier 
transform  of  represents  the  power  spectral  density  (PSD)  of  Sr{t).  Hard¬ 

esty  shows  that  the  shape  of  the  PSD  of  Sr{t)  is  proportional  to  the  probability 
distribution  of  the  radial  wind  velocities  within  the  scattering  volume  (14:175-179). 
Hardesty  and  Menzies  (24:453)  state  that  the  temporal  correlation  time  4  of  the 
received  signal  which  provides  a  measure  of  how  rapidly  Sr{t)  varies  with  time  can 
be  expressed  to  a  rough  approximation  as 

tc^XfcTy  (83) 

where  A  is  the  laser  wavelength  transmitted  and  is  the  standard  deviation  of  the 
distribution  of  aerosol  radial  velocities  within  the  scattering  volume. 
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Chapters  11-13  of  Vaui  Trees’  text  (39:357-536)  provide  additional  insight  into 
the  modeling  required  to  adequately  represent  the  statistical  nature  of  the  backscat- 
tered  signal  Sr{t)  by  using  an  alternate  approach  to  that  presented  in  Equation  (79). 
Van  Trees  provides  a  general,  dense  target  environment  model  which  applies  directly 
to  aerosol  imaging.  Since  the  received  signal  consists  of  a  large  number  of  ran¬ 
domly  distributed  point  target  contributions,  Van  Trees  also  statistically  models  the 
backscattered  signal  as  a  random  process  with  a  time-varying,  Rayleigh  distributed 
amplitude  and  a  time-varying,  uniformly  distributed  phase.  Van  Trees  models  the 
target  backscatter  as  a  complex  Gaussian  random  process  7/3/?  which  can  be  ex¬ 
pressed  in  terms  of  real  and  imaginary  components  as 

niDR  —  'IDRr  +  jlDRi  (84) 

where  '^dRr  and  7£)Hj  are  jointly  Gaussian  random  processes  with  identical  covariance 
functions.  For  the  dense  aerosol  target  case,  the  aerosols  are  assumed  to  compose  a 
continuous  range-spread,  frequency-spread,  speckle  target  with  the  complex  received 
signal  expressed  as 


Sr{t)  =  jEj  /  s{t  -  Td)'rDR{rd,  t  ”  T<i/2)dTa  (85) 

J—oo 


where 

Et  =  Transmitted  energy  (J) 

s{t)  =  Complex  transmitted  signal  (V) 

Sr{i)  =  Sample  function  from  a  complex  Gaussian  random  process  (V) 

Td  =  Range  time-delay  coordinate  with  rj  =  2R/c  (s) 

'iDRijd^i)  =  Complex  Gaussian  random  process  with  independent 
variables  (rd,<)  of  range  and  time. 

Physically,  the  mean  and  covariance  of  'yDRi^di  t  —  Tdl2)  describe  the  statistical  char¬ 
acteristics  of  the  aerosol  target  backscatter.  Van  Trees  makes  two  assumptions  about 
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the  Gaussian  random  process  statistics  of  7£)R(T<j,i  —  T(i/2)  which  are  stated  as  fol¬ 
lows:  Assumption  1.  The  return  from  each  range  interval  Atj,  is  a  sample  function 
of  a  stationary,  zero-mean  complex  Gaussian  random  process.  Assumption  2.  The 
returns  from  different  range  intervals  Arj  are  statistically  independent.  As  a  result 
of  Assumption  1 ,  the  mean  value  o{  ^Dfi{Td,t—Td/2)  \s  equal  to  zero.  And  as  a  result 
of  Assumption  2,  the  covariance  function  of  'fDR{Tdit  —  'rd/2)  designated  as 
can  be  expressed  as 

C^{t^  u)  =  E[yDRiTdi  <)7D/i(‘rd>  w)]  =  KoRiTd,  t  —  u)6{Td  -  r^)  (86) 

where  the  function  KoRi'^^d^i  —  «)  is  dependent  upon  the  reflective  properties  of  the 
target  and  where  the  delta  function  5(Td— t^)  is  a  result  of  the  statistical  independence 
requirement  of  Assumption  2. 

Van  Trees  introduces  the  scattering  function  definition  as 

/OO 

KDR{Td,cx)eyi^{-j2KfdOc)doL  (87) 

•OO 

where  the  scattering  function  SoRird,  fd}  represents  the  Fourier  transform  of  KDRi^d^  <^)j 
therefore,  KDR(Td',oc)  can  be  expressed  as 

/OO 

SDR{rd,  fd)  exp(j2Trfda)dfd  (88) 

•OO 

SDR{Td,fd]  represents  the  range-dependent  PSD  of  the  Gaussian  random  process 
'lDR(Tdi^)-  Since  7uii(Td,t)  is  a  zero-mean,  complex  Gaussian  random  process,  Sr(t) 
as  expressed  in  Equation  (85)  represents  a  sample  function  from  a  zero-mean  complex 
Gaussian  random  process.  The  signal  Sr(i)  at  the  receiver  is  completely  characterized 
by  its  complex  covariance  function  expressed  as 

a,(t,u)  =  E[s,(t)s;(«)]  (89) 
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Substituting  Equation  (85)  into  Equation  (89)  yields  an  expression  for  the  complex 
covariance  function  given  below: 

C'ar(i,  u)  =  EtJ  s{t  -  Tdi)$*{u  ~  Ti)E['fDR{'rd,  t  “  rdl2)'fDR{rd,  t  -  Tj2)]d‘  j 

(90) 

Substitution  of  Equation  (88)  into  Equation  (90)  yields 

/OO 

S{t  -  Td)s’(u  -  T^)KDR{Td,  t  -  u)dTd  (91) 

■OO 

Substitution  of  KDRi^dyCa)  in  terms  of  the  scattering  function  SoRiTdtfd}  given  in 
Equation  (87)  yields  an  alternative  expression  for  Cs^{t^  u)  expressed  as 

u)  =  EtJ  s{t  -  Td)s*{u  -  rd)SDR{rd,  fd]  exp[i27r/d(i  -  u)]dfddTd  (92) 

Recall  that  Sr(t)  represents  the  received  signal  at  the  lidar  receiver  prior  to 
matched-filter  processing.  Equations  (85)  through  (92)  provide  a  summary  of  the 
model  derivation  Van  Trees  (39:357-536)  developed  for  range-spread,  frequency- 
spread  targets  which  is  useful  in  developing  a  scattering  model  for  aerosols.  These 
equations  provide  a  complete  statistical  characterization  of  the  aerosol  backscat- 
tered  signal.  Also,  these  equations  provide  the  foundation  necessary  for  developing 
an  aerosol  target  model  and  obtaining  the  statistical  description  of  the  receiver’s 
MF  response.  As  a  specific  example  of  an  aerosol  scatteiing  function  model  for 
SDR{Tdyfd]y  consider  Shapiro’s  aerosol  target  model. 

Shapiro  (33:1-5)  hcts  developed  an  aerosol  backscatter  model  consistent  with 
the  Van  Trees  dense  target  development  just  described  in  Equations  (85)  through 
(92).  From  the  Shapiro  model,  the  received  signal  from  aerosol  backscatter,  as  in 
Equation  (85),  can  be  expressed  as 

Sr{t)  =  Ej^^  j  S{t  -  Td)‘yDR{Td,  t  -  Tdf2)dTd  (93) 

*/  — OO 
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where 


iDRird,  t)  =  exp(-j27r/df)  /  T{p,  CTd/2,  t  -  CTdl2)d^  (94) 

J^OO 


where 

T{p,Td,t) 

(\p,CTdl2) 

P 

Id 


=  random  speckle  target  field  refiection  coefficient 
=  normalized  spatial  beam  pattern  of  the  transmitter  and 
back  propogated  local  oscillator  at  range  cr<j/2 
=  the  transverse  polar  coordinate  within  the  plane 
perpendicular  to  the  line  of  sight  of  the  radar  beam 
=  the  mean  Doppler  frequency-shift  of  the  aerosols 


Assuming  a  Gaussian  distribution  of  aerosol  relative  velocities,  and  cissuming  that 
the  returns  from  different  spatial  coefficients  p  and  different  range  coefficients  rj  are 
statistically  independent,  the  covariance  of  r(p,  r^,  t)  can  be  expressed  as 


A2 


B[np,u,i)r{-p',r‘,u)]  =  :^^{S^)S(p-p')S(T,-u)pxp 


{t-ul 


(95) 


where 

^{cTdl2)  =  aerosol  backscatter  coefficient  at  range  R  =  CTdl2  {m~^sr~^) 
tc  =  coherence  time  approximated  by  Equation  (83)  (s) 


Substituting  Equations  (94)  and  (95)  into  Equation  (86)  and  solving  for  Kdr{o!^  ^d) 
yields 

\9  y _ VI  y  _ viil  r  ~\ 

(96) 


exp 


Therefore,  the  scattering  function  using  Equation  (87)  can  be  expressed  as 
SDR{rd,  fd]  =  ^{P,  ^)  tc  exp  [-TrHlifd  -  Jdf] 


(97) 
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Van  Trees  (39:361)  defines  the  mean  Doppler  shift  as 

Mdf,  (98) 

where 

2(T^  =  KDRi^di  0)  (99) 

The  scattering  function  5x)fl(Td,  /j)  describes  the  aerosol  range-spread  backscat- 
ter  characteristics  in  terms  of  the  aerosol  backscatter  coefficient  /?(cr<{/2),  and  the 
frequency-spread  backseat  ter  characteristics  in  terms  of  a  Gaussian  distribution  func¬ 
tion.  Note  that  the  smaller  the  coherence  time,  the  broader  the  frequency  spread 
and  therefore  the  broader  the  velocity  distribution  of  the  particles  imaged.  Also 
note  that  the  Gaussian  distribution  term  of  the  scattering  function  is  centered  on 
the  mean  Doppler  frequency-shift  /j.  The  mean  Doppler  frequency-shift  will  vary 
cis  the  pulse  propagates  through  the  atmosphere  and  so  will  the  coherence  time  as 
the  pulse  enters  regions  of  increaising  or  decreasing  turbulence  or  wind  shear.  The 
stationarity  of  the  process  is  limited;  therefore,  this  model  is  valid  over  a  limited 
observation  time. 

When  considering  the  complete  description  of  all  the  point  targets  which  the 
lidar  pulse  encounters  during  its  propogation  through  the  atmosphere,  the  overall 
target  is  seen  to  consist  of  a  continuous  distribution  of  point  target  velocities  and 
ranges.  The  range  can  be  considered  to  extend  beyond  the  maximum  range  of  the 
lidar,  and  the  frequency  will  range  from  zero  to  plus  or  minus  the  highest  Doppler 
frequency-shift  due  to  wind  velocities  encountered.  For  a  mean  radial  wind  velocity 
range  of  0  to  J:40m/s  the  mean  Doppler  frequency-shift  range  is  from  0  to  i:8M Hz. 

In  this  section,  the  signal  Sr{t)  backscattered  from  a  dense  array  of  aerosols 
has  been  modeled  as  a  zero  mean,  complex  Gaussian  random  process.  The  aerosol 
scattering  function  5£)ii{Trf,  fd]  represents  the  range  dependent  PSD  of  the  random 
process.  The  Shapiro  aerosol  model  for  SDR{Td,fd}  represents  a  specific  example  of 
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aerosol  dense  target  modeling  which  is  consistent  with  the  Van  Trees  dense  target 
development.  Radar  measurements  of  SDR{rd,fd}  yield  the  parameters  of  inter¬ 
est  such  as  mean  backscatter  power,  mean  Doppler  frequency-shift,  and  Doppler 
frequency-shift  spread.  Thus  far  only  the  backscattered  signal  Sr{t)  at  the  radar 
receiver  has  been  statistically  modeled.  The  next  logical  step  in  the  development  of 
matched-filter  aerosol  sensing  is  a  statistical  description  of  the  matched-filter  radar 
receiver  response  to  the  backscattered  signal  Sr{t). 


3.5  Conventional  Radar  Receiver 

As  suggested  in  Section  3.2,  the  MF  bank  receiver  response  to  the  signal  Sr{t) 
can  be  expressed  as  the  cross-correlation  with  the  time-translated,  frequency-shifted 
version  of  the  transmitted  signal.  The  discussion  of  resolution  in  Section  3.3  led 
to  the  Rihaczeck  conclusion  that  the  only  means  of  describing  the  MF  response  in 
a  dense  target  environment  of  randomly  distributed  point  targets,  is  through  the 
average  power  response  of  the  matched  filter.  The  average  power  response  was 
shown  to  result  in  a  two-dimensional  convolution  of  the  average  target  cross  section 
in  range  and  Doppler  with  the  signal  ambiguity  function  ^(t,  /).  In  Section  3.4  a 
more  detailed  statistical  description  of  the  aerosol  target  wais  provided.  The  text  of 
Van  Trees  offered  a  general  model  to  describe  the  statistics  of  the  aerosol  target,  and 
a  specific  aerosol  scattering  function  from  Shapiro  was  compared  to  the  Van  Trees 
model.  For  both  the  Rihaczeck  model  and  the  VanTrees  model,  the  mean  power 
response  of  the  MF  represents  the  measureable  quantity  available  to  describe  the 
aerosol  target.  Recall  the  output  of  the  MFSLR  <is  illustrated  in  Figure  3  can  be 
expressed  analytically  as 


o  roo  2 

\y{Tj)\  =  /  Sr{t)s*{t-T)e\p{-j2-Kft)dt 
J—oo 


(100) 


where  Sr{t)  is  now  the  aerosol  backscattered  signal  of  Equation  (85).  The  derivation 
to  determine  the  mean  of  Equation  (100)  is  provided  in  Appendix  A.  This  derivation 
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shows  that  the  mean  MFSLR  response  can  be  expressed  as 

E[\v{rJ)f]  =  N,  +  B,j  j"jDR{rijmr-nJ-U)i!ddri  (101) 

where 

fd  =  target  Doppler  frequency  shift  coordinate  (Hz) 

Td  =  target  range-time  delay  coordinate  (s) 

/  =  MFSLR  frequency  coordinate  (s) 

r  =  MFLSR  range-time  delay  coordinate  (s) 

No  =  mean  Gaussian  spectral  density  (W /Hz) 

The  mean  power  response  of  the  MFSLR  from  Equation  (101)  represents  a  two 
dimensional  convolution  of  the  scattering  function  SDji{Td,fd}  with  the  ambiguity 
function  0{T,f).  Recall  from  Section  3.3  the  comparable  result  for  the  dense  target, 
MFSLR  response  which  Richaczek  models  as  a  two  dimensional  convolution  of  the 
mean  target  cross  section  cijd^fd)  ambiguity  function  0{T^f).  Thus  for 

either  Ccise,  0(t,  /)  along  with  the  target  model,  determine  radar  resolution.  Both 
c{Td,fd)  and  SDR{Td,fd}  represent  PSDs  associated  with  the  backscatter  from  a 
dense  array  of  point  targets  randomly  distributed  in  range  and  velocity.  For  short 
time  duration  pulses  (i.e.  ~  l/is),  where  oijdif^  is  valid,  the  two  models  should 
yield  the  same  analytical  results  in  terms  of  the  MFSLR  response. 

Analytically  modeling  the  aerosol  target  backscatter  with  the  scattering  func¬ 
tion  5£)fl{Td,  fd]  provides  the  advantage  of  statistically  describing  the  aerosol  backscat¬ 
ter  independently  from  the  time  duration  of  the  transmitted  signal.  Thus  the 
derivation  of  the  MFSLR  response  of  Equation  (101)  with  the  scattering  function 
5£)fl{rd,  fd]  represents  a  more  general  solution  where  the  time  duration  of  the  pulse 
can  exceed  the  coherence  time  of  the  aerosol  target.  Radar  resolution  in  terms  of 
Equation  (101)  is  analogous  to  the  radar  resolution  discussed  in  Section  3.3  in  terms 
of  the  Richaczek  dense  target  model.  Thus  one  can  think  of  SuRiTd^fd]  as  a  two 
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dimensional  input,  0(r,  /)  as  the  impulse  response  of  a  two  dimensional  filter,  and 
as  the  two  dimensional  filter  response.  Thus,  the  resolution  with  which 
the  radar  measures  SoRird,  fd]  is  dependent  upon  the  volume  distributions  of  both 
SDR{Tdifd}  and  0(t, /).  This  volume  distribution  cispect  of  radar  resolution  is  dis¬ 
cussed  in  more  detail  in  Chapter  4.  One  can  compare  the  range  zmd  frequency 
resolution  of  a  matched-filter  radar  receiver  to  the  resolution  achieved  with  current 
lidar  systems  employing  Fourier  spectral  analysis. 

As  discussed  in  Section  2.4,  current  lidar  systems  employ  Fourier  spectral  pro¬ 
cessing  techniques  to  obtain  estimates  of  a  range-resolved  power  spectral  density 
called  a  periodogram.  T  samples  of  the  MFSLR  output  are  taken  across  the  fre¬ 
quency  band  for  a  single  range  time  corresponding  to  f  =  t^,  an  estimate  of  the 
PSD  of  the  aerosol  target  is  formed.  The  processing  described  is  comparable  to  the 
current  spectral  processing  technique  discussed  in  Section  2.4.  This  single  estimate 
at  Td  is  range-resolved  according  to  the  resolution  width  At  of  the  ambiguity  func¬ 
tion  0{T,f).  The  lange  of  particles  which  significantly  contribute  to  the  single  PSD 
estimate  at  t  =  tj  for  a  simple  rectangular  pulse  are  contained  within  a  range  swath 
of  width  approximately 

(102) 

This  description  of  range  width  or  resolution  is  comparable  to  the  range  weight¬ 
ing  function  discussed  in  meteorological  radar  literature  (11:52-56).  In  general,  many 
pulses  v/ill  be  noncoherently  integrated  over  the  time  duration  for  which  the  Gaus¬ 
sian  random  process  is  stationary.  Resolution  of  the  PSD  estimate  is  determined  in 
terms  of  the  ambiguity  function’s  frequency  resolution  width  A/.  For  the  simple 
rectangular  pulse.  A/  ~  1/T. 

Note  the  distinction  made  between  resolution  performance  in  terms  of  the 
PSD  estimate  of  SDR{Td,fd}  and  the  resolution  performance  in  terms  of  resolving  a 
single  point  target.  In  general,  single  aerosol  point  targets  are  too  dense  to  resolve 
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individually,  but  if  the  resolution  definition  is  in  terms  of  resolution  with  respect  to  a 
single  point  target,  the  convolution  operation  of  Equation  (101)  on  the  frequency  axis 
must  be  performed.  However  if  resolution  is  defined  in  terms  of  the  resolution  with 
which  SDR{Td,fd}  is  measured,  the  widths  At  and  A/  of  the  ambiguity  function 
will  determine  waveform  resolution  performance.  An  example  will  illustrate  the 
distinction. 

Suppose  as  in  Section  3.3,  that  the  ambiguity  function  can  be  ideally  repre¬ 
sented  as  the  impulse  function 


»(t,/)  =  <(t,/)  (103) 

Substituting  Equation  (103)  into  Equation  (101)  yields  the  expression 

■B[l!/(’-,/)P]=£.SDi!{T,/.}  (104) 

The  form  of  Equation  (104)  provides  the  exact  representation  of  the  scattering 
function  with  infinitesmal  resolution.  However,  if  resolution  is  defined  in  terms  of 
a  single  point  target,  then  resolution  in  frequency  will  be  in  terms  of  the  width  in 
/  of  the  scattering  function  SdrIt,/}  while  resolution  in  range  is  determined  by 
the  width  of  the  ambiguity  function  across  the  t  axis.  This  statement  is  true  as 
long  as  the  MFSLR  response  is  not  coupled  in  frequency  and  range.  If  coupling 
occurs  between  target  range  and  Doppler  (  for  example,  linear  frequency  modulated 
(FM)  pulses),  the  degree  of  frequency  spread  of  the  scattering  function  along  with 
the  ambiguity  function  will  determine  resolution  for  both  frequency  and  range.  In 
the  next  chapter,  signal  waveform  resolution  performance  for  various  types  of  radar 
signals  will  be  determined  in  terms  of  the  scattering  function  SnRiTd^fd}- 
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IV.  Findings  and  Observations 


K 


4.1  Introduction 

As  stated  previously,  the  primary  purpose  of  this  thesis  effort  is  to  compare 
the  resolution  performance  of  a  conventional  matched-filter  radar  receiver  for  various 
signal  waveforms  in  a  dense  target  environment  such  as  aerosols.  It  is  clear  from  the 
derivations  of  Chapter  3  that  resolution  in  a  dense  target  environment  is  determined 
in  terms  of  both  the  scattering  function  and  the  ambiguity  function.  In  general,  the 
backscatter  function  from  Dr.  Shapiro  which  wets  discussed  briefly  in  Section  3.4 
is  assumed  to  adequately  represent  the  aerosol  scattering  function  SDR{Td,fd}-  In 
this  model,  the  target  is  continous  in  range  extending  beyond  the  far  field  range  of 
the  lidar  with  a  reflection  coefficient  determined  by  the  beamwidth  and  the  aerosol 
backscatter  coefficient  j3s{cTdl2).  The  relative  velocities  of  the  aerosols  is  modeled 
in  terms  of  a  Gaussian  distribution  in  frequency  whose  width  is  described  by  the 
coherence  time  U  .  The  coherence  time  in  terms  of  relative  aerosol  velocities  is 
approximately  expressed  as  in  Section  3.3  as  tc  ~  A/ctu.  The  mean  Doppler  frequency 
of  the  aerosols  will  range  from  zero  to  the  highest  expected  atmospheric  velocities. 
Resolution  performance  in  terms  of  the  radar  ambiguity  function  provides  a  measure 
of  the  resolution  of  the  scattering  function  which  the  waveform  can  achieve. 

The  development  of  radar  since  the  early  40’s  has  resulted  in  a  myriad  of 
signal  waveforms  and  signal  processing  techniques,  and  a  detailed  discussion  of  each 
is  well  beyond  the  scope  of  this  paper.  Fortunately  however,  signal  waveforms  may 
be  classified  in  terms  of  their  respective  ambiguity  function  shapes  into  four  distinct 
classes  (31:1078-1083). 

Class  1  waveforms  are  single  pulses  of  carrier  frequency  fc  which  form  an 
ambiguity  function  with  its  volume  concentrated  in  a  central  ridge  at  the  origin  of 
the  (t,/)  plane  with  low  sidelobes.  Class  2  waveforms  also  have  a  concentrated 
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voli  ,vith  low  sidelobes;  however,  the  volume  is  skewed  in  the  (r,  /)  plane.  Single 
pulses  with  linear  frequency  modulation  (FM)  referred  to  as  “chirp”  pulses  represent 
Class  2  waveforms,  and  their  ambiguity  function  is  often  called  the  “knife  edge” 
ambiguity  function.  Class  3  signals  have  an  ambiguity  function  commonly  referred 
to  as  the  “thumbtack”  ambiguity  function.  A  large  variety  of  signal  waveforms 
often  referred  to  as  noise  or  psuedonoise  waveforms  which  require  random  phase  or 
frequency  coding  can  be  clcissified  as  Class  3  waveforms.  In  addition,  single  pulse 
waveforms  with  smoothly  varying  phase  modulation  such  as  quadratic  FM  form  an 
approximate  thumbtack  ambiguity  function.  The  distinctive  characteristic  of  Class 
3  waveforms  is  the  narrow  central  spike  of  negligible  ambiguity  volume  and  the 
broad  pedestal  of  large  B  xT  dimensions  which  contains  most  of  the  volume  of  the 
ambiguity  function.  Class  4  signals  have  an  ambiguity  function  commonly  referred 
to  as  the  “bed  of  nails”.  Periodic  pulse  trains  uniformly  spaced  in  time  represent 
Class  4  waveforms.  Table  1  provides  a  summary  of  the  more  salient  features  of 
each  clciss  of  waveform.  The  four  clcisses  of  waveforms  designated  as  Classes  1,  2, 
3,  and  4  are  discussed  in  Sections  4.2  through  4.5  respectively.  With  these  four 
classifications  of  signal  waveforms  in  terms  of  their  respective  ambiguity  function 
volume  distributions,  the  task  at  hand  is  to  analyze  each  class  of  waveform  with 
respect  to  its  resolution  performance  in  the  aerosol,  dense  target  environment. 


Table  1.  Ambiguity  Function  Parameters  (31:1080) 


Radar  Signals 
Ambiguity  function 

Class  1 

Class  2 

Class  3 

Class  4 

ridge 

sheared  ridge 

thumbtack 

bed  of  nails 

Reolution  cell  size 

unity 

unity 

IITB 

IITB 

Ambiguities 

no 

(r<i,  fd)  coupling 

no 

spikes 

Sidelobes 

low 

low 

high 

low 

TB  product 

unity 

TB>1 

1 

TB  »  1 
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4.2  Class  1  Radar  Signals 


The  single  radar  pulse  of  carrier  frequency  fc  represents  Class  1  radar  wave¬ 
forms.  The  normalized  rectangular  pulse  is  expressed  in  complex  baseband  form 


as 


«<W  = 


(105) 


where  T  represents  the  time  duration  of  the  pulse.  The  ambiguity  function  of  the 
rectangular  pulse  is  derived  from  the  ambiguity  function  definition  of  Section  3.2, 
Equation  (26)  and  is  expressed  as 


x{tJ)  = 


j  (l  -  H)  sine  [(l  -  U)  T/]  if|r|<r 
0  if  \t\  >  T 


(106) 


where  the  function  sinc(a:)  =  (sin7rx)/(7ra;).  If  the  function  tri(T)  is  defined  as 


tri(T)  =  { 


j  (i-M)  it|r|<r 
0  if  |r|  >  T 


(107) 


then  the  ambiguity  function  of  Equation  (106)  can  be  simplified  to 


X(tJ)  =  tri(T)sinc((r/)tri(T)) 


(108) 


Therefore,  the  squared-modulus  ambiguity  function  0(t,  /)  is  expressed  as 


=  tri^(T)sinc^  [(r/)tri(r)] 


(109) 


The  common,  two-dimensional  contour  plot  of  the  ambiguity  function  6{t,  f) 
is  shown  in  Figures  4  and  5  for  long  and  short  time  duration  pulses  respectively. 
The  elliptical  regions  enclose  the  areas  where  the  height  of  ^(t,/)  is  significant 
while  lightly  shaded  regions  depict  the  areas  where  the  height  of  0(t,/)  is  small 
but  nonzero.  The  ellipses  of  Figures  4  and  5  represent  the  approximate  borders 
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of  the  half-power,  equal  height  contours  of  The  ellipse  which  composes 

the  half-power  border  consists  of  all  the  points  of  (t,  /)  which  satisfy  the  equation 
0(r,/)  ~  1/2.  This  convention  for  two-dimensional  contour  plots  introduced  by 
Siebert  (35:217)  simplifies  the  analysis  of  various  waveforms  and  is  common  to  many 
radar  texts  (38,  26,  39).  For  the  rectangular  pulse  ambiguity  function  of  Equation 
(109),  the  shape  of  0(t,/)  along  the  r  and  /  axes  are  expressed  respectively  as 

0(t,O)=  tri^(T)  (110) 

9(0,/)=  sinc’(r/)  (111) 

As  shown  in  Section  3.3,  the  range  time-delay  (RTD)  and  Doppler  frequency- 
shift  (DFS)  resolution  widths  for  the  rectangular  pulse  are  respectively  given  as 

At  ~  r  (112) 

(H3) 

Combining  Equations  (112)  and  (113)  forms  the  nominal  resolution  cell  expressed  cis 

At  A/  ~  1.  (114) 


Similar  results  apply  for  the  Gaussian  pulse  caise.  The  normalized  Gaussian 
pulse  expressed  in  complex  baseband  form  is 


where  T  is  the  root-meaii-square  (RMS)  time  duration  fo  the  signal.  Derivation  of 
the  Gaussian  signal  ambiguity  function  yields 


X('r,/)  =  exp 


f 


(116) 
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Figure  5.  Contour  Plot  of  a  Claiss  1  Ambiguity  Function  (Short  Pulse)  (35:217) 


Therefore,  the  squared-modulus  ambiguity  function  0(r,  /)  can  be  expressed  as 


e(T,/)  =exp 


(117) 


As  in  the  rectangular  pulse  case,  the  equal  height  contours  of  0(t,  /)  form  ellipses, 
and  the  half-power  ellipse  of  the  Gaussian  pulse  is  comparable  to  the  haif-power 
ellipses  of  Figures  4  and  5.  The  shape  of  6{t,  f)  along  the  t  and  /  axes  provides  a 
measure  of  RTD  and  DFS  resolution.  For  the  Gaussian  pulse,  the  shape  of  0(r,  /) 
along  both  axes  is  Gaussian. 

Due  to  the  Fourier  transform  relationship  between  the  time  and  frequency 
domains,  to  a  close  approximation,  the  bandwidth  B  of  any  Class  1  waveform  is 
equal  to  the  reciprocal  of  the  time  duration  T.  Therefore,  Class  1  waveforms  always 
have  nominal  resolution  cells  which  have  the  approximate  dimensions  oi  BxT  with 
an  approximate  area  of  unity.  Figures  4  and  5  illustrate  the  basic  concept.  Figure 
4  represents  the  ambiguity  function  for  a  long  time  duration  pulse,  and  Figure  5 
illustrates  the  ambiguity  function  of  a  short  time  duration  pulse.  As  these  two 
figures  illustrate,  an  increase  in  time  duration  T  results  in  a  decrease  in  the  size 
of  the  frequency  resolution  dimension  A/.  If  the  time  duration  T  is  increrised  to 
infinity,  the  signal  becomes  a  continuous  wave  signal  with  infinitesimal  frequency 
resolution  and  with  no  range  resolution.  Conversely,  if  the  pulse  is  infinitesimally 
small  in  time  duration,  the  impulse  will  provide  infinitesimal  range  resolution  with 
no  frequency  resolution. 

The  coordinate  transformations  to  the  resolution  widths  in  range  and  radial 
velocity  are  respectively  expressed  as 


(118) 

(119) 
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Substitution  of  Equations  (118)  and  (119)  into  Equation  (114)  will  yield  the  nominal 
range-velocity  resolution  cell  given  as 

ARAVr  =  ^  (120) 

For  c  =  3  X  10®  and  A  =  10/im,  ARAV,.  =  750m^/5.  As  Equation  (120)  illustrates, 
the  only  means  of  reducing  the  nominal  range-velocity  resolution  cell  of  a  Class  1 
waveform  is  to  reduce  the  carrier  wavelength.  For  the  simple  Class  1  radar  signal, 
the  resolution  cell  defined  by  Equation  (120)  provides  an  accurate  and  complete  de¬ 
scription  of  the  resolution  with  which  the  radar  measures  the  aierosol  scattering  func¬ 
tion  SDR{Td,fd}-  In  other  words,  the  mean  MFSLR  response  measures  the  aerosol 
scattering  function  with  the  resolution  in  range  and  radial  velocity  of  the  nominal 
resolution  cell  regardless  of  the  distribution  of  SDii{Td,fd}-  The  radar  resolution  of 
the  other  three  classes  of  radar  signal  depend  on  the  volume  distributions  of  both 
SDR{Td,  fd}  and  0(t,  /)  over  the  (r,  /)  plane.  The  remaining  three  classes  of  radar 
signals  to  consider  are  called  pulse  compression  signals.  These  large  time-bandwidth 
product  signals  may  improve  resolution  dependent  upon  the  target  scattering  func¬ 
tion  characteristics.  Recall  from  Section  3.2  that  the  total  volume  of  0(t,  /)  remains 
fixed.  Pulse  compression  radar  signals  simply  redistribute  the  fixed  ambiguity  vol¬ 
ume.  The  first  claiss  of  pulse  compression  signals  designated  as  Class  2  radar  signals 
achieve  large  time-bandwidth  products  through  frequency  modulation. 

4.3  Class  2  Radar  Signals 

Single  pulses  with  linear  frequency  modulation  (FM)  represent  Class  2  signals. 
This  radar  signal  is  commonly  referred  to  in  the  radar  literature  as  a  chirp  pulse,  and 
radars  which  implement  them  are  often  called  chirp  pulse  compression  radars.  The 
technique  of  chirp  pulse  compression  hcis  been  used  extensively  since  its  development 
in  the  40’s,  and  the  literature  describing  its  features  abounds  (38,  26,  39,  30).  As 
Skolnik  explains  (38:420-423),  chirp  pulse  compression  provides  the  radar  designer 
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the  means  of  using  a  long  time  duration  pulse  to  achieve  large  radiated  energy  while 
still  achieving  the  range  resolution  of  a  short  duration  pulse.  In  addition,  means  of 
acquiring  accurate  range  and  velocity  information  from  multiple  pulses  have  been 
developed. 

A  linear  FM  chirp  pulse  can  be  generated  by  either  active  or  paissive  generation 
techniques.  The  active  generation  technique  requires  direct  time-domain  frequency 
modulation  of  the  radar  transmitter.  For  a  rectangular  pulse  of  constant  amplitude 
A  and  time  duration  T  as  illustrated  in  Figure  6a,  the  frequency  modulation  is  linear 
over  the  time  duration  T  of  the  pulse  as  illustrated  in  Figure  6b.  The  slope  \k\  of 
the  linear  FM  can  be  expressed  as 


\k\  =  B/T  (121) 

where  the  bandwidth  B  from  Figure  6b  is  given  as  B  =  /2  —  /i  and  where  T  is  the 
time  duration  of  the  pulse.  The  slope  k  in  general  can  be  positive  or  negative. 

The  pzissive  generation  technique  requires  an  impulse  function  or  short  time 
duration  pulse  as  the  input  to  a  dispersive  phase  filter.  The  frequency  transfer 
function  H{f)  of  the  dispersive  phase  filter  can  be  expressed  in  general  as 

J?(/)  =  l^f(/)lexpy«(/)l  (122) 

where  the  magnitude  of  the  transfer  function  \H{f)\  is  constant  and  the  phase  of  the 
transfer  function  6{f)  is  a  quadratic  function  of  /.  In  general,  the  impulse  response 
of  a  dispersive  filter  with  a  quadratic  phase  function  produces  an  output  signal  whose 
time  domain  phase  is  also  quadratic  and  thus  whose  frequency  varies  linearly  with 
time.  The  response  of  the  dispersive  pheise  filter  to  a  short  time  duration  signal 
with  large  bandwidth  5  is  a  long  time  duration  signal  with  the  bandwidth  B  of  the 
short  pulse.  The  quadratic  phase  of  the  time  domain  signal  output  from  the  filter 
is  equivalent  to  linear  frequency  modulation  as  in  the  active  generation  technique 
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discussed  previously.  Either  the  active  or  the  passive  generation  technique  forms  an 
approximately  equivalent  linear  FM  signal  (30,  38).  For  the  purpose  of  comparing 
the  resolution  performance,  either  generation  technique  will  yield  equivalent  results. 

Consider  the  active  generation  case.  The  transmitted  chirp  pulse  can  be  written 
in  general,  complex  baseband  form  as 

s(t)  =  ls(t)|exp(j7rfcf")  (123) 

where  the  magnitude  |s(t)l  is  a  constant  A  as  shown  in  Figure  6a,  and  the  phase  is 
quadratic  which  is  equivalent  to  the  linear  FM  slope  of  Figure  6b  with  k  =  BfT. 
The  instantaneous  frequency  of  s{t)  is  f{t)  =  kt. 

For  the  analysis  to  follow,  a  detailed  description  of  the  spectrum  of  the  signal 
s{t)  of  Equation  (123)  is  unnecessary.  The  Fourier  transform  of  the  general  chirp 
pulse  of  Equation  (123)  can  be  written  in  general  form  as 

S(/)  =  |5(/)Nxp  [;•«{/))  (124) 

where  0{f)  is  a  quadratic  function  of  frequency  and  the  shape  of  |5(/)|  is  a  function 
of  both  the  envelope  of  the  transmitted  signal  and  the  bandwidth  B  representing  the 
degree  of  linear  frequency  modulation  (30,  3).  If  the  time-bandwidth  product  TB 
of  the  transmitted  signal  is  large  {TB  1),  it  can  be  shown  that  the  bandwidth 
of  the  spectrum  envelope  is  approximately  the  linear  FM  modulation  bandwidth  B 
where  B  =  f 2  —  fi  as  shown  in  Figure  6b  (38:423).  Often  B  is  called  the  •:hirp 
bandwidth.  Thus  the  transmitted  chirp  pulse  has  an  envelope  of  time  duration  T 
and  an  increased  bandwidth  of  approximately  B. 

To  compress  the  pulse  in  the  time  domain,  a  matched  filter  is  employed.  The 
frequency  transfer  function  can  be  expressed  in  general  cis 

H(.f)  =  SV)  =  |S(/)I  exp  l-m)]  (125) 
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Amplitude 


Figure  6.  Amplitude  (a)  and  Frequency  (b)  of  the  Transmitted  Chirp  Pulse  (38:423) 
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where  the  transfer  function  H{f)  is  the  complex  conjugate  of  the  transmitted  signal 
spectrum  written  in  general  form  in  Equation  (124).  Thus  the  matched-filter  is 
represented  as  a  dispersive  phase  filter  with  the  quadratic  phase  characteristic  —0{f). 
As  a  result,  the  filter  removes  the  ’  onlinear  phase  6{f)  of  the  received  signal  and 
compresses  the  received  signal  into  a  short  duration  pulse.  The  degree  of  compression 
of  the  MF  response  is  proportional  to  the  bandwidth  B.  Although  a  large  time- 
bandwidth  signal  is  transmitted,  the  MF  response  compresses  the  time  duration  to 
form  a  pulse  of  duration  IjB  with  a  bandwidth  of  jB;  thus,  the  MF  response  produces 
a  signal  with  a  time-bandwidth  product  of  unity  {TB  Oi  1)  (30:55).  Derivation  of  the 
chirp  pulse  ambiguity  function  will  provide  a  more  precise  description  of  the  chirp 
pulse  MF  response. 

Consider  the  active  generation  of  a  rectangular-envelope,  chirp  pulse.  The 
transmitted  signal  can  be  written  in  normalized,  complex  biiseband  form  cis 

(126) 

Similarly,  the  Gaussian-envelope  transmitted  signal  can  be  written  in  normalized 
complex  b2iseband  form  as 

where  T  is  the  RMS  signal  time  duration.  The  rectangular  and  Gaussian  enve¬ 
lope  pulses  represent  the  two  most  common  pulses  shapes  considered  in  radar  texts. 
The  Gaussian  chirp  pulse  yields  a  slightly  broader  MF  response  without  significant 
sidelobes  eis  compared  to  the  rectangular  chirp  pulse.  For  both  these  pulses,  the 
linear  FM  ambiguity  functions  are  similar  in  form  and  interpretation,  so  only  the 
rectangular  pulse  ambiguity  function  is  analyzed  in  detail. 
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One  of  the  properties  of  the  ambiguity  function  (39:290)  can  be  stated  as 
follows:  If 

s{t)<r^x{T'J)  (128) 

where  the  symbol  means  “has  the  ambiguity  function”,  then 

s(<)  exp  i-*  xC'T,  /  -  kr)  (129) 

Therefore,  the  ambiguity  function  of  the  chirp  pulse  is  simply  the  ambiguity  function 
of  the  unmodulated  pulse  of  the  same  time  duration  with  its  major  axis  rotated  to 
the  slope  of  f  =  kr  as  illustrated  with  a  two-dimensional  contour  plot  of  6{t,  f)  in 
Figure  7.  Substituting  the  expression  for  the  ambiguity  function  of  the  unmodulated 
rectangular  pulse  in  Equation  (108)  into  Equation  (129),  the  chirp  pulse  ambiguity 
function  can  be  written  as 

X('r,  /)  =  tri(r)sinc  (T(/  -  A:T)tri(r)]  (130) 

where  again  \k\  =  B/T.  From  the  definition  of  the  sine  function.  Equation  (130)  can 
be  expressed  as 

sin  [7rr(/  -  fcr)tri(T)) 

XKT,J)  TrT{f-kT) 

For  large  time-bandwidth  product  signals  {TB  ^  1),  the  tri(T)  term  in  the  nu¬ 
merator  of  Equation  (131)  represents  a  slowly  varying  distortion  term  which  can  be 
ignored  (30:171).  Therefore,  the  approximate  envelope  of  the  ambiguity  function  is 
expressed  as 

X(r,  /)  ~  sine  [ttT  (/  -  kr)]  (132) 

Thus  the  squared-modulus  ambiguity  function  0{r,f)  is  given  in  approximate  form 
as 

0(r,  /)  ~  sinc^  —  /:t)]  (133) 
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Figure  7,  Contour  Plot  of  a  Class  2  Ambiguity  Function  (35:217) 
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Equations  (132)  and  (133)  illustrate  the  symmetry  of  the  ambiguity  function  in  r 
and  in  /.  A  vertical  slice  parallel  to  either  axis  yields  a  sine  function  MF  response. 
This  symmetry  indicates  the  coupling  effect  between  target  range  and  velocity  for 
chirp  pulse  compression  radar.  Recall  from  Section  3.2  that  the  ambiguity  function 
represents  the  response  of  an  array  of  MFs  shifted  in  frequency.  An  alternative  view 
of  the  ambiguity  function  is  also  possible.  The  ambiguity  function  also  represents 
the  response  of  a  single  MF  as  the  Doppler  frequency  of  the  received  signal  is  varied 
(38:411-412).  Due  to  the  symmetry  of  Equations  (132)  and  (133),  the  MF  response 
for  a  Doppler  mismatch  fd  is  equivalent  to  the  MF  response  translated  by  Atd  =  fdjk 
for  fd  =  0.  As  a  result,  a  MF  receiver  for  the  linear  FM  chirp  pulse  employs 
a  single  MF.  This  technique  is  possible  since  Doppler  frequency  shifts  within  the 
broad  bandwidth  B  of  the  filter  simply  produce  an  additional  time  lead  or  time  lag 
(At^  =  fdjk)  on  the  t  axis  dependent  upon  the  slope  of  k  employed  with  limited 
amplitude  attenuation  or  distortion  (30:248). 

For  the  single  MF  case,  letting  /  =  0  in  Equation  (133)  yields  the  MF  response 
on  the  T  axis  expressed  as 

0(t,  0)  =  sinc^  (ttT^t)  (134) 

Substituting  k  =  BfT  into  Equation  (134)  yields 

0(r,  0)  =  sinc^  (ttBt)  (135) 

Therefore  as  discussed  previously,  the  single  MF  response  of  Equation  (135)  com¬ 
presses  the  received  \«fnal  to  a  time  duration  of  1/B  with  the  broad  bandwidth  B. 
As  a  result,  the  large  ne-bandwidth  product  of  the  received  signal  is  compressed 
in  time  to  a  MF  output  signal  with  T5  ~  1, 

In  summary,  the  MF  receiver  for  linear  FM  chirp  pulses  consists  of  a  lineai 
FM  frequency  characteristic  which  is  opposite  to  the  frequency  characteristic  of  the 
transmitted  chirp  pulse.  The  receiver  compresses  the  received  signal  from  a  time 
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duration  T  to  a  time  duration  1/jB  with  a  bandwidth  of  B.  The  MF  output  is  a 
compressed  signal  with  TB  c:i  1.  The  chirp  pulse  receiver  can  employ  a  single  MF 
since  a  Doppler  frequency  shift  of  the  received  signal  within  the  bandwidth  B  of 
the  filter  produces  an  additional  time  lead  or  time  lag  in  the  MF  response  given  as 
Atj  =  fd/k  which  is  dependent  upon  the  slope  of  k  of  the  transmitted  signal  with 
limited  attenuation  or  distortion. 

The  time-frequency  coupling  of  the  linear  FM  signal  allows  a  considerable 
simplification  of  the  radar  receiver  since  only  one  MF  is  required.  However,  since 
both  the  target’s  range  and  the  target’s  radial  velocity  cause  delays  in  the  MF 
response,  the  receiver’s  simplicity  also  results  in  a  limitation  of  the  receiver’s  ability 
to  determine  both  target  range  and  radial  velocity  simultaneously  from  one  pulse. 
In  fact,  this  receiver  achieves  no  resolution  in  Doppler  frequency  shift  within  the 
bandwidth  B. 

Before  resolution  properties  of  the  linear  FM  signal  are  discussed,  first  consider 
the  simple  case  where  the  radar  images  a  single  point  target  of  unknown  RTD 
and  unknown  DFS  fd.  As  a  result  of  the  target  DFS  /j,  the  output  peak  Tp  of  the 
MF  response  occurs  at 

Tp  =  Td  +  ^Td  (136) 

where  Atj  =  fd/k.  Thus,  for  a  moving  target,  the  peak  of  the  MF  response  does 
not  represent  the  true  target  range  tj.  The  symmetry  of  Equations  (132)  and  (133) 
shows  that  all  point  targets  with  combinations  of  tj  and  fd  which  lie  along  the  ridge 
of  the  ambiguity  function  satisfy  the  following  equation: 

fd  -kTd  =  0  (137) 

All  combinations  of  {Td,fd)  which  satisfy  Equation  (137)  will  produce  a  peak  at  the 
same  location  Tp.  Therefore,  the  range  measurement  Tp  is  ambiguous  over  the  time 
duration  T  of  the  chirp  pulse.  Klauder  (21:804)  and  Rihaczek  (30:246-248)  show  that 


64 


there  is  a  technique  called  “range  extrapolation”  which  resolves  the  range  ambiguity 
inherent  to  a  chirp  pulse  compression  receiver  imaging  a  single  point  target.  As 
Rihaczek  explains,  if  the  point  target  DFS  fd  is  converted  to  the  radial  velocity 
coordinate  VJ.,  the  range  error  due  to  At^  can  be  expressed  as 


AR  = 


(138) 


Since  the  target  is  moving  with  the  radial  velocity  VJ.,  the  target  travels  a  distance 
of  AR  in  At  =  ARfVr  seconds.  Therefore,  the  range  error  AR  of  Equation  (138) 
results  in  a  travel  time  At  expressed  as 


(139) 


Rihaczek  states  that  as  a  result  of  Equation  (139)  the  range  indicated  by  the  delay 
of  the  MF  response  peak  corresponds  to  the  range  not  at  the  time  of  signal  impact 
but  rather  the  range  at  time  At  later  or  earlier  depending  on  the  direction  of  the 
linear  frequency  sweep  k.  Therefore,  the  linear  FM  receiver  provides  a  measure  of 
extrapolated  range  Rg.  The  extrapolated  range  Re  for  a  positive  slope  k  is  measured 
At  after  signal  impact  and  Rg  for  a  negative  slope  k  is  measured  At  before  signal 
impact  (30:247-248).  The  measurement  of  extrapolated  range  Re  provides  the  means 
of  accurately  determining  target  range;  otherwise,  unambiguous  range  measurements 
with  a  chirp  pulse  compression  receiver  is  limited  to  stationary  targets. 

The  chirp  pulse  compression  receiver  can  accurately  extrapolate  the  range 
through  the  simple  time  translation  At  =  fc/k  which  is  a  constant  of  the  radar  even 
though  the  target  DFS  fd  is  unknown.  However,  the  receiver  provides  no  means  of 
accurately  extrapolating  the  target  velocity  from  a  single  pulse.  When  both  target 
RTD  and  DFS  are  unknown,  the  chirp  pulse  compression  receiver  cannot  simultane¬ 
ously  determine  both  the  target  paramters  and  fd  from  a  single  pulse.  The  target 
DFS  fd  can  be  determined  only  from  consecutive  extrapolated  range  measurements 
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(30:249).  Converting  to  range  and  radial  velocity  coordinates  and  assuming  the  tar¬ 
get  velocity  does  not  significantly  change  between  consecutive  range  measurements, 
the  extrapolated  target  velocity  can  be  expressed  as 


(140) 


where  ile,  and  Re^  represent  the  consective  extrapolated  range  measurements  from 
two  independently  processed  chirp  pulses  seperated  in  time  by  Tp.  At  this  point, 
the  relevant  issues  of  measuring  the  range  and  radial  velocity  of  a  single  point  target 
with  a  chirp  pulse  compression  receiver  have  been  discussed.  This  discussion  provides 
an  introduction  to  the  topic  of  chirp  pulse  radar  resolution  which  is  the  next  topic 
presented. 

For  the  case  where  numerous  point  targets  are  simultaneously  imaged  rather 
than  a  single  point  target,  the  resolution  capability  of  the  chirp  pulse  radar  becomes 
the  primary  concern.  Since  the  chirp  pulse  radar  measures  only  extrapolated  range 
with  no  capability  of  measuring  radial  velocity  from  a  single  pulse,  point  targets  with 
the  same  extrapolated  range  cannot  be  resolved  (30:251-252).  The  discussion  with 
respect  to  Equations  (136)  and  (137)  showed  that  all  point  targets  with  combinations 
{Td,fd)  which  satisfy  Equation  (137)  and  thus  lie  along  the  central  ridge  of  the  chirp 
pulse  ambiguity  function  will  produce  a  MF  peak  at  the  same  location  Tp  on  the  r 
axis.  As  a  result,  when  multiple  point  targets  have  the  RTD  and  DFS  combinations 
which  satisfy  Equation  (137),  all  the  point  target  responses  of  the  MF  output  sum 
together  at  r  =  Tp  to  form  a  single  peak  MF  response  at  Tp  which  is  inseperable. 
Therefore,  in  a  dense  target  environment  where  the  point  targets  are  spread  in  both 
range  and  radial  velocity,  the  degree  of  range  resolution  which  the  chirp  pulse  radar 
can  achieve  is  a  function  of  both  the  width  of  the  MF  response  which  is  approximately 
1/B  from  Equation  (135)  and  the  width  of  the  range  and  velocity  “spread”  of  point 
targets  which  cause  them  to  lie  along  the  central  ridge  of  the  chirp  pulse  ambiguity 
function. 
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Consider  a  scenario  where  good  range  resolution  is  the  primary  requirement  of 
the  radar.  Since  range  resolution  performance  of  the  chirp  pulse  radar  is  dependent 
upon  the  distribution  in  range  and  radial  velocity  of  the  point  targets  which  compose 
the  dense  target  environment,  one  can  consider  three  cases  of  dense  target.  In  each 
case  the  dense  target  is  cleissified  in  terms  of  the  velocity  behavior  of  the  point 
targets.  For  the  first  case,  the  point  targets  are  stationary.  For  the  second  case,  all 
the  point  targets  move  with  the  same  radial  velocity  And  finally  for  the  third 
case,  the  point  targets  move  with  a  mean  radial  velocity  Vr  and  move  relative  to 
each  other.  Each  dense  target  case  is  discussed  in  the  following  paragraphs. 

The  first  czise  is  the  easiest  to  analyze.  Since  all  the  point  targets  are  stationary, 
no  range-Doppler  coupling  occurs.  For  this  cjise  the  scattering  function  SDR{Td,fd} 
can  be  expressed  as 

SDR{Td,fd}  ==  SDR{Td]S{fd)  (141) 

Since  only  a  single  matched  filter  is  employed  for  linear  FM,  the  frequency  variable 
/  can  be  set  to  zero.  Thus,  the  mean  MFLSR  response  from  Section  3.5  can  be 
written  as 

E  [|j/  (r,  0)f]  =  No +  EtJ  Sdr  {t^,  fd]  0  (r  -  tj,  -fd)  dfddrd  (142) 

Substituting  Equation  (141)  into  Equation  (142)  yields 

E  [|y  (r,-0)p]  =  No  ^Et  Sdr  {tj}  0  (r  -  r^,  0)  dvd  (143) 

The  form  of  Equation  (143)  shows  that  no  coupling  exists  between  the  target  range 
and  Doppler  frequency.  Therefore,  the  dense  stationary  target  is  simply  a  range- 
spread  target  where  the  width  of  0{t  —  Td,  0)  across  the  r  axis  represents  the  width 
of  the  MF  response  in  RTD  At.  The  RTD  resolution  width  At  which  the  chirp 
pulse  MFSLR  can  achieve  is  approximately  1/B,  and  the  range  resolution  width  AR 
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is  expressed  as 


(144) 

Thus,  the  width  of  6[t  —  rj,  0)  across  the  r  axis  determines  the  range  resolution  of 
the  chirp  pulse  radar. 

Next,  consider  the  second  case  where  all  the  point  targets  move  at  the  radial 
velocity  VJ.  and  are  stationary  with  respect  to  each  other.  The  scattering  function 
fd}  can  be  expressed  as 

5'Dn{T<i,  fd)  =  SDR{rd}S{fd  —  fdr)  (145) 

where  fd^  is  the  DFS  of  the  point  targets  due  to  the  radial  velocity  VJ..  Substitution  of 
Equation  (145)  into  the  expression  for  the  mean  MFSLR  response  given  in  Equation 
(142)  yields 

E  [ijr  (t.  /)!"]  =  JV.  +  SoR  {Ti,  U]  9  (r  -  TJ,  -/*)  dn  (146) 

As  in  the  first  czise,  the  form  of  Equation  (146)  shows  that  no  coupling  exists  between 
the  target  range  and  Doppler  frequency.  If  range  extrapolation  is  employed,  the 
width  of  6(t,  f)  across  the  r  axis  determines  the  range  resolution  width  At.  As  in 
Case  1,  the  RTD  resolution  width  At  is  approximately  1/B,  and  the  range  resolution 
A/2  is  expressed  as 

(147). 

In  terms  of  velocity  measurement,  the  method  of  velocity  extrapolation  as  given  in 
Equation  (140)  will  not  work  here  because  the  range  is  continuous,  so  the  radar  has 
no  means  of  resolving  with  respect  to  velocity.  Thus,  the  receiver  provides  no  means 
of  extrapolating  the  radial  velocity  Vy  because  the  target  is  continuous  in  range. 

And  finally  consider  the  third  case  where  the  point  targets  move  with  a  mean 
velocity  Vy,  and  point  target  relative  motion  is  nonzero.  For  the  third  case,  range- 
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Doppler  coupling  in  occurs.  This  third  case  represents  the  aerosol  target  en¬ 
vironment  which  is  consistent  with  the  Shapiro  aerosol  scattering  model  given  in 
Section  3.4.  Due  to  the  range-Doppler  coupling  associated  with  chirp  signals  and 
illustrated  by  the  chirp  ambiguity  function  of  Figure  7,  the  degree  of  point  target 
relative  motion  will  determine  the  degree  of  range  resoluton  which  the  chirp  radar 
can  achieve.  A  specific  example  in  terms  of  Dr  Shapiro’s  aerosol  target  model  will 
provide  a  more  quantitative  description  of  the  range  resolution  At  which  the  chirp 
radar  can  achieve. 

Dr.  Shapiro  has  demonstrated  the  range  resolution  potential  of  a  chirp  pulse 
radar  for  aerosol  imaging  (33:1-5).  Recall  from  Section  3.4  the  Shapiro  model  of  the 
aerosol  target  and  from  Section  3.5  the  mean  power  MFSLR  response.  The  mean 
power  MFSLR  response  was  expressed  in  general  cis 


E  [Ij;  (t,  /)|“]  =  ^  J  >•(*)«■(<  -  ’■)  exp  * 


(148) 


where  r{t)  =  Sr{t)  -1-  n{t)  and  Sr{t)  was  expressed  as 


Sr{t)  =  [  s{t  -  Td)'rDR{Td,  t  -  Tdl2)dTd  (149) 

J— 00 


and  7Dii(7'<i,t)  which  represents  the  Shapiro  aerosol  model  was  expressed  as 

7Dfl(Td,  t)  =  exp{-j2irfdt)  I  T{p,,  CTd/2,  t  -  Td/2)  {p,  CTd/2)  dp  (150) 

*/— 00 


To  simplify  the  analysis,  Dr  Shapiro  assumes  that  the  mean  Doppler  frequency- 
shift  (DFS)  fd  is  equal  to  zero.  Since  the  range  can  be  extrapolated  in  the  more 
general  case  where  the  mean  DFS  is  nonzero,  the  simplification  of  setting  the  mean 
DFS  to  zero  is  not  especially  limiting.  With  this  simplification,  a  single  matched 
filter  with  /  =  0  is  employed.  Thus  the  mean  power  MFSLR  response  of  Equation 
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(148)  simplifies  to 


■  B[ls/(T,0)fl=£ 


r{t)s*{t  —  T)dt 


|2l 


(151) 


Substituting  Equation  (149)  and  Equation  (150)  into  Equation  (151)  yields  an  an¬ 
alytical  expression  for  the  mean  power  MFLSR  response.  The  detailed  calculations 
to  determine  the  analytical  expression  of  the  mean  power  MFSLR  response  to  an 
aerosol  target  are  laborious  and  unnecessary  for  the  present  discussion  of  range  res¬ 
olution.  Dr  Shapiro’s  results  are  given  without  proof;  however,  Dr  Shapiro’s  results 
were  verified  using  the  general  model  for  the  mean  power  MFSLR  response  with 
/  =  0.  Therefore,  the  mean  MFSLR  response  of  Equation  (142)  can  be  expressed  as 


E  [|y  (2" >  0)P]  =  No +  Et  J  j  ^  Sdr  {Td,  fd}  0{t-  Td,  —fd)  dfddtd  (152) 


Recall  that  Sbr  {Td,  fd)  represents  the  PSD  of  the  complex  Gaussian  random  process 
7r>R(Td,t)  of  Equation  (150)  where  /d  =  0.  Also  recall  from  Section  3.4  the  Shapiro 
aerosol  model  where  Sdr  {rd,  fd]  is  expressed  as 


Sdr  {Td,  fd]  =  ^  [f>,  y)  (-^^^cfd)  (153) 


Substituting  Equation  (153)  into  Equation  (152)  yields  a  general  result  for  the  mean 
power  MFSLR  response  in  terms  of  the  Shapiro  aerosol  target  model  and  the  ambi¬ 
guity  function  of  the  transmitted  radar  signal.  Again,  the  detailed  calculations  are 
unnecessary.  Either  Equation  (152)  or  Equation  (151)  yield  the  same  expression  for 
the  mean  MFSLR  resp^.'-  ^  the  aerosol  target  model.  The  results  are  stated  as 
follows: 
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Dr  Shapiro  showed  that  transmission  of  a  Gaussian  chirp  pulse  expressed  in 
complex  baseband  form  as 

[-  *"]  (^54) 

will  result  in  a  mean  power  MFSLR  response  expressed  in  general  form  as 

E  [ly(r,  of]  =  Q(Td)VF(r  -  Td)dTi  (155) 


where  Q(Td)  is  the  infinite  resolution,  average  aerosol  backscatter  profile,  and  where 
W{t  —  Td)  is  a  Gaussian  range  wieghting  funcvion  which  determines  the  range  reso¬ 
lution  of  the  chirp  pulse  radar.  W{t  —  Td)  can  be  expressed  in  general  form  as 


W{t  —  Td)  =  K  exp 


8(t  -  Tdf 
At 


(156) 


where  K  is  a  system  constant,  and  At  is  the  full  width  in  the  r  coordinate  between 
the  attenuation  points  of  iy(r  —  Td)-  Shapiro  shows  that  Ar  can  be  expressed 
as 


At  = 


2l/2y 


(157) 


where  T  is  the  RMS  time  duration  and  B  is  the  RMS  bandwidth  of  the  transmitted 
signal  s{t).  Therefore,  At  represents  the  RTD  resolution  of  the  chirp  pulse  radar. 
Multiplying  Equation  (157)  by  c/2  yields  an  analogous  expression  for  the  range 
resoluton  AR  given  as 


AR  = 


_ cTj2^l^ _ 

[l  +  (irBr/4)V(l+3’V«?)]'''‘ 


(158) 


For  more  details  of  the  derivation,  refer  to  Dr  Shapiro’s  report  (33:1-5). 

The  primary  interest  in  terms  of  range  resolution  of  chirp  pulse  radar  is  the 
nature  of  the  range  resolution  AR  of  Equation  (158).  Recall  that  the  coherence 
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time  tc  is  a  measure  of  the  relative  velocities  of  the  aerosol  particles.  Therefore, 
as  stated  earlier,  the  degree  of  range  resolution  AR  which  the  chirp  pulse  radar 
achieves  is  a  function  of  the  distribution  of  cierosol  radial  velocities  the  lidar  beam 
illuminates.  When  the  point  targets  such  ais  aersols  move  relative  to  each  other, 
the  scattering  function  is  broadened  in  frequency.  As  a  result,  the  range  resolution 
AR  is  broadened  to  the  extent  of  significant  overlap  of  the  chirp  pulse  ambiguity 
function  0(t,/)  and  the  scattering  function  SDR{Td^fd]  as  illustrated  in  Figure  8. 
The  schematic  of  Figure  8  is  consistent  with  the  results  of  Equations  (158).  Both 
show  that  AR  is  a  function  of  the  volume  distributions  of  0(t,/)  and  SDR{Td,fd}> 


Figure  8.  Combined  Contour  Plots  of  the  Scattering  and  Ambiguity  Functions 


As  an  example,  let  T  =  5/is,  U  =  1/xs,  and  B  =  lOMHz.  Substitution  into 
Equation  (158)  will  yield  AR  ~  72.6m.  If  the  coherence  time  were  extremely  large 
which  indicates  minimal  relative  motion  of  aerosols,  AR  ~  27m.  Also  in  comparison 
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to  an  unmodulated  pulse  of  the  same  time  duration,  AR  ~  cr/2  ~  750m.  It  is  clear 
from  this  example  that  with  the  assumption  of  a  fixed  mean  velocity  of  aerosols,  the 
relative  velocity  distribution  of  the  aerosols  somewhat  degrades  the  pulse  compres¬ 
sion  of  the  matched-filter  but  still  improves  range  resolution  performance  compared 
to  an  unmodulated  pulse  of  the  same  time  duration.  To  achieve  the  same  range  res¬ 
olution  with  the  unmodulated  pulse,  the  pulse  duration  T  required  is  approximately 
equal  to  0.5  microseconds.  Due  to  transmitter  peak  power  limitations,  a  5  /is  pulse 
may  reduce  the  maximum  range  of  lidar  operation  .  As  in  the  second  dense  target 
Ccise,  no  means  of  extrapolating  the  mean  velocity  is  available  since  the  target  is 
continuous  in  range. 

In  summary,  three  cases  of  dense  target  have  been  considered,  and  for  each 
case,  a  single  linear  FM  matched-filter  was  employed.  For  Cases  1  and  2  where  the 
point  targets  were  assumed  stationary  with  respect  to  each  other,  the  width  of  0(r,  0) 
across  the  r  axis  determined  the  range  resolution  AR  of  the  chirp  pulse  radar.  For 
Case  3,  assuming  =  0  to  simplify  the  analysis,  relative  motion  of  the  point  targets 
Wcis  included.  As  a  result,  the  range  resolution  width  AR  became  a  function  of  the 
volume  distributions  of  both  the  scattering  function  Sqr  {tj,  fd]  and  the  chirp  pulse 
ambiguity  function  0(t,/).  For  Cases  2  and  3  in  general,  the  chirp  pulse  MFSLR 
could  not  extrapolate  the  mean  velocity  because  the  aerosol  target  was  continuous 
in  range. 

To  cover  the  topic  of  chirp  pulse  compression  radar  as  thoroughly  as  possible, 
an  additional  approach  for  extrapolating  target  range  and  velocity  of  chirp  pulses  is 
presented.  This  second  approach  of  chirp  radar  processing  '"ported  by  Ramp  and 
Wingrove  (29:75-82)  requires  the  transmission  of  two  pulses  of  alternate  linear  FM 
slopes  and  corresponding  matched-filters.  Each  signal  is  processed  independently 
and  combined  to  provide  an  interpretation  of  point  target  range  and  velocity.  The 
basic  signal  processing  for  a  single  point  target  can  be  explained  as  follows:  The 
transmission  of  a  chirp  pulse  with  positive  slope  k  will  result  in  a  compressed  pulse 
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MF  response  which  peaks  at 


T(+)  =  Td-  ^Td  (159) 

where  Atj  =  fdf  |fc|.  Similarly,  for  a  chirp  pulse  of  negative  slope  A:,  the  MF  response 
peaks  at 

T(_)  =  Td  +  Ard  (160) 

Therefore,  the  true  target  range-time  delay  Td  can  be  found  by  adding  Equations 
(159)  and  (160)  to  obtain 

_  T(+)  +  T(-) 

2 

Also,  subtracting  Equation  (159)  from  Equation  (160)  and  substituting  Arj  =  fdf  |A:| 
yields  the  expression  for  the  DFS  fd  given  as 

A  =  (162) 

Thus  Equation  (161)  measures  the  target  RTD  and  Equation  (162)  measures 
the  target  DFS.  This  type  of  signal  processing  described  by  Ramp  and  Wingrove  has 
been  demonstrated  with  a  CO2  laser  radar  (16:152-162).  While  this  type  of  receiver 
signal  processing  accurately  determines  the  range  and  the  velocity  of  a  single  point 
target  without  having  to  rely  upon  range  extrapolation,  this  processing  technique 
can  determine  neither  range  nor  radial  velocity  of  point  targets  in  a  dense  target 
environment  such  as  aerosols.  As  Rihaczek  explains,  with  multiple  point  targets  it 
will  generally  be  impossible  to  order  the  received  returns  into  their  respective  pairs 
to  determine  their  combined  range  and  radial  velocity  (30:251). 

It  has  been  shown  that  pulse  compression  radar  will  provide  improved  range 
resolution  perforrnance  in  terms  of  an  aerosol  target  as  a  function  of  the  degree  of 
relative  velocity  of  the  aerosols  at  the  expense  of  obtaining  no  velocity  informa¬ 
tion.  If  the  radar  designer  desires  both  range  and  velocity  information  of  a  target 
which  consists  of  point  targets  distributed  in  range  and  velocity,  the  chirp  pulse 


does  not  offer  a  viable  alternative  to  the  simple  'inmodulated  pulses  designat  *  ’  .  ? 
Class  1  radar  signals.  The  remaining  two  classes  »■.'  pulse  compression  radar  signals 
to  consider  require  less  detailed  analysis.  The  next  case  to  consider  is  the  large 
time-bandwidth  product,  pulse  compression  waveform  which  forms  the  thumbtack 
cimbiguity  function. 

4.4  Class  3  Radar  Signals 

A  variety  of  signal  waveforms  often  referred  to  as  noise  or  pseudonoise  wave¬ 
forms  with  an  ambiguity  function  shape  which  approximates  i’i3  ambiguity  diagram 
of  Figure  9  are  designated  as  Class  3  radar  signals  (38:418).  This  idealized  ambiguity 
function  is  commonly  referred  to  as  the  thumbtack  ambiguity  function  in  obvious 
reference  to  its  shape.  The  approximate  thumbtack  ambiguity  function  of  Figure  9 
consists  of  a  narrow  central  spike  containing  a  negligibly  small  portion  of  the  am¬ 
biguity  volume.  The  nominal  resolution  cell  is  defined  in  terms  of  the  central  spike 
disjoint  from  the  base  of  the  ambiguity  function  and  is  expressed  as 

(163) 

where  T  is  the  RMS  time  duration  of  the  transmitted  signal,  and  B  is  the  RMS 
bandwidth  of  the  tra*ismitted  signal.  The  volume  of  unity  distributed  across  the 
area  of  dimensions  TxB  with  an  approximate  height  oiljTB  is  called  the  pedestal 
of  the  ambiguity  function.  Large  time-bandwidth  product  TB'>\  waveforms  which 
consist  of  iong  time-duration,  large  bandwidth  signals  generate  the  thunibtack  am¬ 
biguity  function.  An  infinite  variety  of  signals  referred  to  as  noise  waveforms  such  as 
signals  obtained  through  random  phase  or  frequency  coding  generate  approximate 
thumbtack  ambiguity  functions.  Another  signal  waveform  which  provides  a  poorer 
approximation  to  the  thumbtack  ambiguity  function  is  waveforms  with  smoothly 
varying  frequency  modulation  such  as  quadratic  FM  (31:«')'  1  ]07./j  More  precise, 
three-dimensional  ambiguity  diagrams  for  various  Clzus  3  waveforms  n.  >  be  found 


ArA/  = 


TB 


75 


in  the  radar  literature  (30). 

As  Rihaczek  (30:133-140)  and  Siebert  (35:219-220)  explain,  the  narrowness  of 
the  central  spike  of  the  thumbtack  ambiguity  function  provides  excellent  range  and 
velocity  resolution  for  the  siijgle  target  as  long  as  the  signal  can  be  detected  above 
the  self  clutter  power  interference.  The  self-clutter  is  obtained  from  contributions 
due  to  the  point  targets  within  the  broad  pedestal  of  dimensions  Tx  B,  If  the  cross 
section  of  the  target  within  the  central  spike  of  the  ambiguity  function  is  significantly 
stronger  than  the  cumulative  cross  section  of  all  the  point  targets  within  the  pedestal 
of  the  ambiguity  function,  then  the  Class  3  waveform  will  provide  improved  resolution 
performance.  To  quantify  this  concept,  radar  texts  usually  introduce  the  term  S/C 
for  signal  to  clutter  ratio. 

Recall  from  Section  3.5  that  the  mean  power  MFSLR  response  is  a  two- 
dinmensional  convolution  of  '  he  scattering  function  Sdr  {t^,  fd}  with  the  ambiguity 
function  5(r,  /)  which  in  this  case  is  the  thumbtack  ambiguity  function.  The  bound¬ 
aries  where  SoRifd^fd]  is  significant  (i.e.  e“^  attenuation  width)  in  the  (t,/)  plane 
defines  the  occupied  target  space.  Since  the  scattering  function  is  continuous  in 
range,  the  range  extent  of  the  scattering  function  always  extends  beyond  the  dimen¬ 
sion  T  of  the  pedestal  across  the  r  axis.  If  /j  is  designated  eis  the  boundary  where 
Sdr  fd}  is  significant  across  the  /  axis,  then  either  fb  <B  ox  fh>  B  where  B  is 
the  dimension  of  the  pedestal  across  the  /  axis.  The  analysis  of  the  two  cases  where 
either  fh>  B  ox  fb  <  B  are  briefly  summarized  from  the  Rihaczek  text  (30:133-140). 

First  consider  the  Ccise  where  fb  >  B,  The  desired  nominal  resolution  cell  ie 
defined  in  terms  of  the  widths  across  the  t  and  /  axes  of  the  ambiguity  function’s 
central  spike  as  in  Equation  (163).  The  pedestal  S'mface  of  the  ambiguity  function 
can  be  divided,  by  the  area  of  the  nominal  resolution  cell  to  determine  the  number 
N  of  resolution  cells  within  the  pedestal.  Thus,  the  pedestal  contains  N  nominal 
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resolution  cells  where  N  is  expressed  as 


TV- 


TB 

1/TB 


~  {TBf 


(164) 


Let  a  represent  the  value  of  the  scattering  function  within  a  resolution  cell.  If  (Xc 
represents  the  cumulative  cross  section  of  the  scattering  function  from  each  resolu¬ 
tion  cell,  then  the  mean  cross  section  of  the  scattering  function  in  each  cell  can  be 
expressed  approximately  as 


— 


N 


(165) 


Therefore,  if  cr,  represents  the  cross  section  value  of  the  scattering  function 
Sdr  {Td,  fd}  at  the  resolution  cell  within  the  central  spike  of  the  ambiguity  func¬ 
tion,  the  average  S/C  can  be  expressed  as 


S  <x,{TB) 

G-  (JmN 


(166) 


Substituting  for  N  from  Equation  (164)  into  Equation  (166)  yields 

c  \tbJ  KvJ 


(167) 


If  the  minimum  signal  to  clutter  ratio  S/C  is  arbitrarily  chosen  as  unity,  then  Equa¬ 
tion  (167)  constrains  the  ratio  cxs/o'm  to 


—  c^TB  (168) 

Om. 

Equation  (168)  shows  that  to  detect  the  signal  with  a  minimum  signal  to  clutter 
ratio  of  unity,  a,  must  be  TB  greater  than  average.  In  other  words,  to  achieve  the 
resolution  of  the  nominal  resolution  cell  defined  in  Equation  (163)  with  S/G  —  1,  the 
scattering  function  cross  section  a,  must  be  TB  higher  than  the  average  cross  section 
level.  Thus,  any  resolution  cell  of  the  scattering  function  for  which  <  {TB  X  Cm) 
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is  completely  masked  in  the  self-clutter. 

Next,  consider  the  second  case  where  fb  <  B.  For  this  case,  resolution  perfor¬ 
mance  will  improve  somewhat.  It  is  easily  shown  that  when  fb  <  5,  Equation  (167) 
can  be  expressed  as 


As  in  the  first  case,  for  a  signal  to  clutter  ratio  of  unity.  Equation  (169)  cam 
be  expressed  ais 

—  ^  Tfb  (170) 

The  form  and  interpretation  of  Equation  (170)  is  analogous  to  that  of  Equation 
(168).  To  achieve  the  resolution  of  the  nominal  resolution  cell  defined  in  Equation 
(163)  with  S/C  =  1,  the  scattering  function  cross  section  cr,  must  be  Tfb  greater 
than  average.  Thus  the  resolution  improves  dependent  upon  the  degree  of  significant 
overlap  between  fb  and  B. 

For  either  case,  it  is  clear  that  the  thumbtack  ambiguity  function  does  not 
achieve  the  desired  resolution  of  the  nominal  resolution  cell.  In  general,  portions  of 
the  scattering  function  where  <7,  is  approximately  less  than  or  equal  to  the  mean 
value  a-jn  are  completely  masked  by  the  self-clutter;  therefore,  in  a  dense  target 
environment,  the  true  resolution  which  Class  3  signals  achieve  is  defined  by  the  broad 
dimensions  of  the  ambiguity  function  pedestal.  Rihaczek  (30:133-140)  shows  that  in 
general,  the  resolution  of  Class  3  signals  in  a  dense  target  environment  is  inferior  to 
the  resolution  of  a  simple  unmodulated  pulse  designated  as  a  Cl2iss  1  waveform.  For 
these  reasons,  Rihaczek  and  Siebert  state  that  the  inferior  resolution  performance 
of  Class  3  radar  signals  make  them  in  general  a  poor  choice  of  waveform  in  a  dense 
target  environment  (30;  35).  Therefore,  Class  3  radar  signals  are  not  a  viable  option 
in  comparison  to  Class  1  radar  signals.  The  final  radar  signal  to  consider  forms  the 
so-called  bed  of  nails  ambiguity  function. 
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Ji.S  Class  4  Radar  Signals 

Class  4  radar  signals  form  the  so-called  “bed  of  nails”  ambiguity  function. 
A  uniformly-spaced  coherent  pulse  train  as  in  Figure  10  illustrates  the  common 
example  of  a  Class  4  waveform.  The  periodic  pulse  train  is  essentially  a  means  of 
producing  a  long  time-duration  signal  for  improved  frequency  resolution  with  a  short 
time-duration  pulse  which  provides  good  range  resolution.  If  p(i)  represents  a  single 
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Figure  10.  Periodic  Pulse  Train 


complex  pulse  of  duration  T  with  the  ambiguity  function  designated  as  Xpi'’’}  f)  >  ^he 
complex  periodic  pulse  train,  normalized  with  respect  to  its  energy  can  be  expressed 

cLS 

1  M-l 

=  JHJl  E  P(^  “ 

771=0 

where  Tr  is  the  pulse  repetition  period  and  M  represents  the  total  number  of  pulses. 
To  siiiiplify  the  ana:lysis,  the  pulse  train  can  be  centered  at  time  t  =  0  so  that  p{t) 
is  expressed  as 

■* 

1  iW-l  ^  1  \ 

4t)  =  +  iiM  -  l)r.)  (172) 

The  pulse  train  of  Equation  (172)  is  illustrated  in  Figure  10  for  the  common  pulse 
duration  of  M  =  3.  Using  Equation  (172)  to  determine  the  ambiguity  function, 
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Blahut  (8:6-17-6-18)  has  shown  that  the  ambiguity  function  can  be  expressed  eis 


X{rj)  = 


1  Af^l 

¥  S 

Tn=-(A/-1) 


s)n(7r/j[r) 


If  it  is  assumed  that  the  pulses  are  of  short  duration  in  comparison  to  the  pulse 
repetition  interved  JV,  the  components  of  the  ambiguity  function  Xpi'^^  f)  which  are 
centered  at  r  =  mTr  will  overlap  only  with  their  low  sidelobe  parts,  so  to  a  close 
approximation,  one  can  look  at  each  component  ambiguity  function  individually. 
The  component  ambiguity  function  for  one  specific  value  of  m  can  be  written  as 


Xm(r,/) 


M  —  |m| 
M 


Xp{T-mTrJ) 


sin7r/(M  —  ImDTr 
{M  —  \m\)smTrfTr  _ 


(174) 


The  first  term  of  Equation  (174)  within  brackets  forms  a  broad  triangular  envelope 
across  the  t  axis.  Within  this  triangular  envelope,  the  ambiguity  function  XpiT^f) 
is  periodically  repeated  at  each  value  of  m  across  the  r  axis.  Across  the  /  axis,  Xp 
forms  a  broad  envelope,  and  the  term 


sin7r/(M  —  |m|)rr 
sin  lifTr 


(175) 


samples  the  broad  frequency  envelope  of  Xp  f  —  ^/Tr  where  k  is  an  integer.  The 
end  result  is  an  ambiguity  function  cis  illustrated  in  Figure  11  with  fine  structure 
both  across  the  t  axis  and  across  the  /  axis  (30:187). 

To  determine  the  nominal  resolution  cell  dimensions  in  terms  of  the  ambigu¬ 
ity  function’s  central  spike,  the  widths  across  each  axis  at  the  origin  are  measured 
according  to  any  of  the  three  resolution  criteria.  For  RTD  resolution,  let  m  =  0  and 
f  —  0  to  simplify  Equation  (174)  to 


Xm=o(r,0)  =  xp(r,0) 


(176) 
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Figure  11.  Contour  Plot  of  a  Class  4  Ambiguity  Function  (35:217) 
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which  represents  the  ambiguity  function  across  the  r  axis  of  a  single  pulse  of  duration 
T;  therefore,  the  RTD  resolution  width  At  ~  T  ~  1/R  is  the  same  as  for  a  single 
pulse  p{t)  of  duration  T  of  the  coherent  pulse  train.  Similarly,  to  determine  frequency 
resolution  of  the  central  spike  across  the  /  axis,  let  m  =  0  and  t  =  0  so  that 


Xm=o(0)  /)  —  /) 


sin  TT  fMTr 
sin  irfTr 


(177) 


The  first  null  of  Equation  (177)  occurs  at  /  =  l/MTr  =  l/Tj;  therefore,  the  DFS 
resolution  width  can  be  written  as 


A/  = 


MTr  Td 

The  nominal  time-frequency  resolution  cell  is 

1 


(178) 


At  A/  a 


BTd 


(179) 


Range  and  radial  velocity  resolution  widths  can  be  expressed  respectively  cis 


AR  ~ 

C 

w 

(180) 

AVr  ^ 

A 

2Td 

(181) 

And  therefore,  the  nominal  range- velocity  resolution  cell  can  be  expressed  cis 


cA  1 


ARAVr  ~  —  X 


TdB 


(183) 


Thus,  the  coherent  pulse  train  provides  the  range  resolution  of  a  single  pulse 
of  duration  T  and  the  improved  velocity  resolution  of  a  single  pulse  of  duration 
Td-  The  cost  of  this  improved  resolution  is  range  and  velocity  ambiguities  as  il¬ 
lustrated  in  Figure  11  where  again  the  ambiguity  function  ^(t,/)  is  approximated 
with  a  two-dimensional  contour  plot.  Range  ambiguities  are  seen  to  occur  at  multi- 
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pies  of  the  pulse  repetition  interval  Tr  and  frequency  ambiguities  are  seen  to  occur 
at  multiples  of  l/Tr.  These  range  and  frequency  ambiguities  are  approximately  of 
the  same  magnitude  as  the  central  spike;  therefore,  if  targets  are  present  beyond 
Tr  in  range  or  l/Tr  in  frequency,  significant  contributions  from  targets  outside  the 
nominal  resolution  cell  will  degrade  resolution  in  either  domain.  Typically,  coher¬ 
ent  pulse  trains  are  employed  when  the  range  and  velocity  extent  of  the  target  or 
targets  of  interest  is  confined  to  the  immediate  region  surrounding  the  central  spike 
of  the  ambiguity  function  as  illustrated  in  Figure  12  (30:143).  When  the  targets 
are  confined  to  this  region,  the  MFSLR  response  is  of  no  interest  beyond  the  value 
Ti  corresponding  to  the  outer  limit  of  target  range  and  /i  for  target  velocity.  To 
create  a  large  “clear”  rectangular  area  about  the  central  spike,  one  would  logically 
desire  that  both  range  and  velocity  ambiguities  are  seperated  as  much  as  possible. 
However,  since  the  spread  in  range  ambiguity  is  determined  by  Tr  and  the  spread  in 
velocity  ambiguity  is  determined  by  l/Tr,  a  compromise  between  range  ambiguity 
and  velocity  ambiguity  is  obviously  required.  The  maximum  unambiguous  range  is 

^  (183) 

and  the  maximum  unambiguous  velocity  is 

Kn  =  (184) 

The  product  of  the  maximum  unambiguous  range  and  maximum  unambiguous  ve¬ 
locity  is  a  constant  of  the  radar  expressed  as 

=  ±J  (185) 

Equation  (185)  illustrates  the  trade-off  required  between  maximum  unambiguous 
range  and  maximum  unambiguous  velocity.  This  limitation  is  especially  troublsome 
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for  the  atmospheric  sensing  lidar  operating  at  a  wavelength  of  10/xm.  Substituting 
A  =  10{im  and  c  =  3  x  10®m/s  gives 

w? 

iCKn  =  ±750—  (186) 

s 

Equation  (186)  shows  that  the  product  BmVm  is  simply  too  small  for  most  mete¬ 
orological  targets  of  interest  using  a  CO2  lidar.  For  example,  the  CTI  lidar  with 
an  aerosol  range  of  approximately  5  kilometers  will  consist  of  velocity  ambiguities 
beyond  0.0375m/s.  Or  conversely,  the  CTI  lidar  discussed  in  Section  2.4  can  un¬ 
ambiguously  detect  aerosols  of  velocities  from  zero  to  ±35m/s.  For  14i  =  ±35m/s, 
Rm  —  21.4m  These  results  show  that  for  the  case  of  aerosol  imaging,  coherent  pulse 
trains  are  not  in  general  a  viable  option  for  improving  combined  range  and  velocity 
resolution. 

In  summary,  each  of  the  four  classes  of  radar  signals  have  been  evaluated 
with  respect  to  their  resolution  capability.  The  ambiguity  function  for  each  claiss 
of  waveform  was  evaluated  in  terms  of  the  resolution  with  which  the  target  scat¬ 
tering  function  SoRi^d,  fd]  can  be  determined.  The  pulse  compression  waveforms 
designated  as  Classes  2,  3,  and  4  were  each  compared  to  the  resolution  capability 
of  the  simple  unmodulated  pulse  designated  its  a  Cleiss  1  radar  signal.  In  each  case, 
the  pulse  compression  waveforms  did  not  provide  improved  resolution  performance 
in  both  range  and  radial  velocity  simultaneously.  As  a  result,  none  of  the  pulse 
compression  waveforms  evaluated  represent  a  viable  alternative  to  the  Class  1  radar 
signal  currently  in  use  in  lidar  systems. 


V.  Conclusions  and  Recommendations 


Results  of  this  thesis  effort  provide  a  thorough  review  of  conventional  matched- 
filter  radar  theory  as  it  applies  to  an  aerosol  sensing  lidar.  Chapter  2  provided  an 
historical  perspective  of  lidar  development  and  a  brief  discussion  of  current  lidar 
systems  and  signal  processing  techniques.  Chapter  3  provided  a  discussion  of  the 
baisic  theory  of  conventional  matched-filter  radars,  the  complex  aerosol  target  model, 
and  the  MFSLR  response  to  the  aerosol  target  model.  The  development  of  the 
MFSJ'R  response  given  in  Appendix  A  is  considered  to  be  an  original  result  which 
provider  ‘he  MFSLR  response  to  any  complex  target  satisfying  the  conditions  of  the 
Van  Trees  target  model  given  in  Section  3.4.  In  Chapter  4,  each  of  the  four  classes  of 
radar  signals  were  evaluated  in  terms  of  their  resolution  performance.  Each  class  of 
radar  signal  was  evaluated  in  terms  of  the  resolution  with  which  it  could  reproduce 
the  scattering  function  of  the  target. 

5.1  Conclusions 

As  a  result  of  this  thesis  effort,  the  mean  MFSLR  response  to  a  dense,  aerosol 
target  environment  can  be  expressed  as 

^  (|2/(t, /)|]  =  No +  EtJ Sonird,  fd}e (r  -  r^,  /  -  /<,)  drddfd  (187) 

where  SDR{Td)fd}  represents  the  scattering  function  of  the  aerosol  target  or  any 
complex  target  which  satisfies  the  conditions  of  the  Van  Trees  target  model  given  in 
Section  3.4  and  where  0(t— tj,  / — fd)  represents  the  squared-modulus,  matched-filter 
response  to  a  single  point  target  with  range  time-delay  Td  and  Doppler  frequency- 
shift  fd.  The  MFSLR  response  derived  in  Appendix  A  represents  a  two-dimensional 
convolution  of  the  target  scattering  function  Sonird,  fd}  with  the  signal  ambiguity 
function  6{T,f).  The  resolution  with  which  the  MFSLR  measures  the  scattering 
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function  is  dependent  upon  the  volume  distributions  of  both  the  ambiguity  function 
and  the  aerosol  scattering  function.  Since  the  normalized  volume  of  0(t,  /)  is  always 
fixed  to  equal  unity,  the  four  classes  of  radar  waveforms  are  represented  by  ambiguity 
functions  with  different  distributions  of  the  same  volume  over  the  (t,  /)  plane.  In 
each  case,  the  waveform’s  resolution  capability  with  respect  to  the  target  scattering 
function  was  evaluated.  The  pulse  compression  waveforms  designated  as  Class  2,  3, 
and  4  radar  signals  did  not  provide  an  improved  resolution  capability  in  comparison 
to  simple  unmodulated  pulses  designated  as  Class  1  radar  signals  currently  in  use 
with  aerosol  sensing  lidars.  In  conclusion,  if  a  conventional  MFSLR  is  employed. 
Class  1  radar  signals  represent  the  best  choice  in  terms  of  resolution  performance 
and  receiver  design  simplicity. 

5.2  Recommendations 

The  question  of  whether  more  sophisticated  pulse  compression  radar  signals 
provide  iinproved  resolution  performance  with  a  conventional,  matched-filter  radar 
heis  been  thoroughly  researched,  and  the  coi.  elusion  is  that  they  do  not.  Class 
1  radar  signals  represent  the  best  choice  signal  for  aerosol  sensing  lidars  when  a 
conventional,  matched-filter  receiver  is  employed.  However,  this  thesis  effort  did  not 
determine  whether  the  MFSLR  represents  the  optimum  receiver  design  for  a  dense, 
aerosol  target  environment.  It  may  be  possible  to  make  better  use  of  the  statistics 
of  the  aerosol  target  with  a  more  sophisticated  receiver  design.  The  many  receiver 
designs  which  Dr  Van  Trees,  develops  in  his  text  (39)  may  offer  some  insight  into 
more  sophisticated  receiver  designs  to  achieve  better  resolution  performance. 
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Appendix  A.  Derivation  of  the  MFSLR  Response 


Consider  the  schematic  of  Figure  3  which  illustrates  an  array  of  matched-filters 
followed  by  square-law  envelope  detectors.  The  mean  power  response  of  the  MFSLR 
can  be  expressed  analytically  cis 


R  [l2/(7’>  f)?]  -R  [^>-(0  +  ^(01 -  t)  exp  {-j27ift) 


dt 


(188) 


Equation  (188)  can  be  written  as 


R[\y{r,f)f]  =  j  J_^E{[sr{t)  +  n{t)][s;{u)  +  n*{u)]} 

X  [s’*  {i  —  t)s{u  —  t)]  exp  [—j2‘K f  {t  —  u)]  dtdu 


=  jj  E[sr(t)s*(u)-f  Sr(t)n*(u)-l-s*(ti)n(i)-|-n(<)n’'(u)] 
xs*(i  —  t)s  {u  ~  t)  exp  [— i27r f  {t  —  w)]  dtdu  (189) 

where  n{t)  if.  a  zero-mean  complex  white  Gaussian  noise  process  due  to  the  shot 
noise  limited  photodetector  whose  covariance  function  Cn(t,u)  is  expressed  as 

Cn{t,u)  =  E[n(t)n*(M)] 

=  No6{t  -  u)  (190) 

where  No  is  the  spectral  density  hu/v  of  the  photodetector  shot  noise  limited,  white 
Gaussian  noise  process.  Taking  the  expected  value  of  each  term  within  brackets  of 
Equation  (189)  ,  the  cross- terms  E  [sr(0^*(^)]  R  [^?(^)^(0]  S®  since  n{t) 

is  a  zero-mean  random  process.  Therefore,  Equation  (189)  can  be  written  as 
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=  j  J  ^E[sr{i)s;{u)]s*{i-T)s{u-T)exT;)[-j2i:f{t-u)]didu 
+  J  J  E[n{t)n*{u)]s^{t-T)s{u-T)exp[-j27rf{t-u)]dtdu 

(191) 

Substituting  Equation  (190)  into  Equation  (191)  and  simplifying  yields 

■E[b(’'./)P]  =^‘'  +  /  J  ^C,^t,u)s'{t  -  t)s(u  -  r)exp[-i2-xf  {t  -  u)]dtdu 

(192) 

where  Csr  is  the  complex  covariance  of  the  received  signal  expressed  as 

Cs,{i,u)  =  E[sr{i)s;{u)]  (193) 

and  where  Sr{t)  has  been  normalized.  Recall  from  the  Van  Trees  model,  the  expres¬ 
sion  for  the  received  signal  given  again  eis 

Sr(t)  =  E^^  I  s{t  -  Td)'i/DRird,t  -  ~)dTd  (194) 

where  'yDR{Td,i)  is  a  zero-mean  complex  Gaussian  random  process  whose  complex 
covariance  C^{i^u)  is  expressed  as 

C^{i,u)  =  E[rfDR{Td,t)7hR{'^d^'^)] 

=  KDR{jd,t-u)S{Td-Td)  (195) 

Substitution  of  Equation  (194)  into  Equation  (193)  yields 

r  roo  f  f 

Cs,{t,u)  =  EiJ  C^{t,u)s{t  -  Td)s*{u  -  Td)dTddT^  (196) 
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Substituting  Equation  (195)  into  Equation  (196)  and  simplifying  yields 

/oo 

—  Td)s*{u  -  Td)KDR(Td,  t  -  u)dTd  (197) 

-00 

Recall  from  Section  3.4,  the  Fourier  transform  relation  between  KoRiTdtOi)  of  Equa¬ 
tion  (197)  and  the  scattering  function  SDR{Tdifd}  expressed  as 


/<x> 

SoRiTd,  fd}  exp(j27r/da)d/d 

-OO 


(198) 


Substituting  Equation  (198)  into  Equation  (197)  yields 

^sr{t,u)  =  s{t  -  rd)s*(u  -  Td)SDR  {rdyfd} 

X  exp  [j2Trfd{t  -  ii)]  dTddfd  (1: 

Substituting  Equation  (199)  into  Equation  (192)  yields 

E  [l2/('r,  f)f]  =  /  /^  -  fd){i  -  «)] 

X  /  /  s{t-Td)s’'{u-Td)SDR{Tdifd}dTddfddtdu 
J  J^CO 


(200) 


Rearranging  terms,  Equation  (200)  can  be  expressed  as 


K  +  j"jDR{u,h} 

/OO 

$"{1  -  T)s{t  -  Td)  exp  [-i27r(/  -  fd)t]  di 

-OO 

/CO 

s*(w  —  t)s(«  —  Td)  exp  [— i27r(/  —  fd)u]  dudrddfd 

-OO 


=  No  +  EtjJ  SDR{Td,fd} 

X  [  s{t)$*[t-{T-Td)]exp[-j2ir{f-fd)t]< 

J— OO 
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Recall  the  variable  translations  of  Section  3.2  which  led  to  the  ambiguity  funct’ 
definition.  The  variable  translations  were 


T  =  T  -  Td 


f  =  f-fd 


(202) 

(203) 


Substituting  Equations  (202)  ana  (203)  into  Equation  (201)  yields 


N(,  +  Et  j  I  Sdr  {Tdi  fd} 

/°°  I  I 

s(t)s*(t  —  r  )  exp{—j2'irf  t)dt 

■CO 

f<X>  ,  , 

X  /  s’‘{u)s{u  —  T  )exp{j2'iTf  u)dudTddfd  (204) 


Recall  the  ambiguity  function  definition  given  as 

x{tJ)=  [  s{t)s*{t-T)exTp{-j2Trft)dt  (205) 

J  —  OO 

where  the  prime  notation  has  been  dropped  for  convenience.  Thus,  the  squared 
modulus  of  x('^,  f)  designated  as  0(t,  f)  can  be  expressed  as 


=  (  s(t)s*{t  —  t)  Qxp{—j27rft)dt 

J  —  OO 
rco 

X  /  s{u)s*{u-  t)cxp(  j2T:fu)du  (206) 

J-OO 


Substituting  Equation  (206)  into  Equation  (204)  yields  the  final  expression  for  the 
MFLSR  response  given  as 

E[\y{rJ)f]=No  +  Et  j  J"^SDR{TdJd]  6  {r  -  r^J  -  fd)dT,dfd  (207) 

Therefore,  the  mean  power  response  of  the  MFSLR  is  in  terms  of  the  pho¬ 
todetector  shot  noise  energy  No  and  a  two-dimensional  convolution  of  the  scattering 
function  SDR{Td,fd}  with  the  radar  signal  ambiguity  function  0(t, /).  Equation 
(207)  is  useful  for  evaluating  the  MFSLR  response  in  general  for  an  arbitrary  radar 
ambiguity  function  and  arbitrary  distributed  target  whose  backscattered  signal  is 
described  as  a  zero-mean,  complex  Gaussian  random  process. 


93 


Bibliography 


1.  Bachman,  C.  G.  Laser  Radar  Systems  and  Techniques.  Dedham,  Mass:  Artech 
House,  1979. 

2.  Barnes,  Norman  P.  and  Charles  E.  Byvik.  “Mid-Infrared  Solid-State  Lasers 
and  Laser  Materials,”  Proceedings  of  SPIE  Airborne  and  Spaceborne  Lasers  for 
Terrestrial  Geophysical  Sensing,  889:  92-107.  Bellingham,  WA;  SPIE,  1988. 

3.  Barton,  David  K.  Radars,  Volume  3.  Dedham  MA:  Artech  House,  1975. 

4.  Biernson,  G.  and  R.  F.  Lucy.  “Requirements  of  a  Coherent  Laser  Pulse-Doppler 
Radar,”  Proceedings  of  the  IEEE,  51:  202-213  (January  1963). 

5.  Bilbro,  J.  and  others.  “Airborne  Doppler  Lidar  Wind  Field  Measurements,” 
Bulletin  American  Meterorological  Society,  65:  348-359  (April  1984). 

6.  Bilbro,  James  W.  “Atmospheric  laser  Doppler  velocimetry:  an  overview,”  Op¬ 
tical  Engineering,  19:  533-542  (July  1980). 

7.  Bilbro,  James  W.  and  William  W.  Vaughn.  “Wind  Field  Measurement  in  the 
Nonprecipitous  Regions  Surrounding  Severe  Storms  by  an  Airborne  Pulsed 
Doppler  Lidar  System,”  Bulletin  American  Meteorological  Society,  59:  1095- 
1100  (September  1978). 

8.  Blahut,  Richard  E.  Theory  of  Remote  Surveillance.  Reading  MA:  Adison- 
Wesley,  1988. 

9.  Byer,  Robert  L.  “Diode  Laser-Pumped  Solid-State  Lasers,”  Science,  239:  742- 
747  (February  1988). 

10.  Curran,  Robert  J.  “Satellite-Borne  Lidar  Observations  of  the  Earth:  Require¬ 
ments  and  Anticipated  Capabilities,”  Proceedings  of  the  IEEE,  IT.  478-490 
(March  1989). 

11.  Doviak,  Richard  J.  and  Dusan  S.  Zirnic.  Doppler  Radar  and  Weather  Observa¬ 
tions.  Orlando  FL:  Academic  Press,  1984. 

12.  Goodman,  Joseph  W.  Statistical  Optics  New  York:  John  Wiley  and  Sons,  1985. 

13.  Hardesty,  R.  Michael  and  others.  “A  ground-based,  injection-locked,  pulsed 
TEA  laser  for  atmospheric  wind  meeisurements,”  Proceedings  of  SPIE  Airborne 
and  Spaceborne  Lasers  for  Terrestrial  Geophysical  Sensing,  889:  23-28.  Belling¬ 
ham,  WA:  SPIE,  1988. 

14.  Hardesty,  Michael  R.  Measurement  of  Range-Resolved  Water  Vapor  Concentra¬ 
tion  by  Coherent  CO2  Differential  Absorption  Lidar.  PhD  dissertation.  Naval 
Postgraduate  School,  Monterey  CA,  1984  (AD-A143  594). 

15.  Hardesty,  Michael  R.,  Senior  Research  Scientist.  Telephone  interview.  NOAA 
Wave  Propogation  Laboratory,  Boulder  CO,  12  August  1990. 


94 


16.  Henderson,  D.  M.  and  others.  “Pulse  Compression  of  an  FM  Chirped  CO2 
Lciser,”  Proceedings  of  SPIE  Fifth  Conference  on  Coherent  Laser  Radar:  Tech¬ 
nology  and  Applications,  1181:  66-78.  Bellingham,  WA;  SPIE,  1989. 

17.  Huffaker,  R.  M.  and  others.  “Laser-Doppler  System  for  Detection  of  Aircraft 
Trailing  Vortices,”  Proceedings  of  the  IEEE,  58:  322-326  (March  1970). 

18.  Huffaker,  R.  Milton.  “Performance  Analysis  and  Technical  Assessment  of  Co¬ 
herent  Lidar  Systems  for  Airborne  Wind  Shear  Detection,”  Proceedings  of  SPIE 
Airborne  and  Spaceborne  Lasers  for  Terrestrial  Geophysical  Sensing,  889:  65- 
76.  Bellingham,  WA:  SPIE,  1988. 

19.  Kavaya,  Michael  J.  and  others.  “Remote  wind  profiling  with  a  solid-state 
Nd:YAG  coherent  lidar  system,”  Optics  Letters,  If-  776-778  (August  1989). 

20.  Kavaya,  Michael  J.  and  Robert  T.  Menzies.  “Lidar  aerosol  backscatter  mea¬ 
surements:  systematic,  modeling,  and  calibration  error  considerations,”  Applied 
Optics,  Sf:  3444-3453  (November  1985). 

21.  Klauder,  J.  R.  “The  Theory  and  Design  of  Chirp  Radars,”  The  Bell  System 
Technical  Journal,  Vol.  XXXIX  No.  4:  745-808  (April  I960). 

22.  Lawrence,  I.  P.  and  R.  M.  Huffaker.  “A  Lciser  Velocimeter  for  Remote  Wind 
Sensing,”  The  Review  of  Scientific  Instruments,  fS:  512-518  (March  1972). 

23.  Lawrence,  T.  R.  and  others.  “Feasibility  and  design  considerations  of  a  global 
wind  sensing  coherent  infrared  radar  (WINDSAT),”  Proceedings  of  SPIE 
Physics  and  Technology  of  Coherent  Infrared  Radar,  300:  34-43.  Bellingham, 
WA:  SPIE,  1981. 

24.  Menzies,  Robert  T.  and  R.  Micheal  Hardesty.  “Coherent  Doppler  Lidar  for 
Measurements  of  Wind  Fields,”  Proceedings  of  the  IEEE,  77:  449-462  (March 
1989). 

25.  Menzies,  Robert  T.,  Senior  Research  Scientist.  Telephone  interview.  Caltech  Jet 
Propulsion  Laboratory,  Pasadena,  CA,  10  July  1990. 

26.  Nathanson,  Fred  E.  Radar  Design  Principles.  New  York:  McGraw-Hill,  1969. 

27.  Petheram,  J.  F.  and  others.  “The  Leiser  Atmospheric  Wind  Sounder- 
Preliminary  Design,”  Proceedings  of  SPIE  Fifth  Conference  on  Coherent  Laser 
Radar:  Technology  and  Applications,  1181:  66-78.  Bellingham,  WA:  SPIE, 
1989. 

28.  Post,  M.  J.  and  others.  “National  Oceanic  and  A.tmospheric  A.dministration’s 
(NOAA)  pulsed,  coherent,  infrared  Doppler  lidar  -  characteristics  and  data,” 
Proceedings  of  SPIE  Fifth  Conference  on  Coherent  Laser  Radar:  Technology 
and  Applications,  300:  60-65.  Bellingham,  WA:  SPIE,  1989. 

29.  Ramp,  H.  0.  and  E.  R.  Wingrove.  “Principles  of  Pulse  Compression,”  IRE 
Transactions  on  Military  Electronics,  MIL-5:  109-116  (April  1961). 


95 


30.  Rihaczek,  August  W.  Principles  of  High-Resolution  Radar.  New  York:  McGraw- 
Hill,  1969. 

31.  Rihaczek,  August  W.  “Radar  Waveform  Selection  -  A  Simplified  Approach,” 
IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  AES-7,  No.  6:  1078- 
1086  (November  1971). 

32.  Schwiesow,  R.  L.  and  R.  E.  Cupp.  “Remote  Doppler  velocity  measurements  of 
atmospheric  dust  devil  vortices,”  .Applied  Optics,  15:  1-2  (January  1976). 

33.  Shapiro,  Jeffrey  H.  Coherent  Laser  Radar  Remote  Sensing,  Appendix  B,  Fi¬ 
nal  Report,  19  January  1984.  Massachusetts  Institute  of  Technology,  Lincoln 
Laboratory. 

34.  Siebert,  William  McC.  “A  Radar  Detection  Philosophy,”  IRE  Transactions  on 
Information  Theory,  IT-2,  No.  3:  204-221  (September  1956). 

35.  Siebert,  William  McC.  “A  Radar  Detection  Philosophy,”  IRE  Transactions  on 
Information  The..-  y,  TT-2,  No.  3:  204-221  (September  1956). 

36.  Sirmans,  Dale  and  Bill  Bumgarner.  “Numerical  Comparison  of  Five  Mean  Fre¬ 
quency  Estimators,”  Journal  of  Applied  Meteorology,  14:  991-1003  (September 
1975). 

37.  Skolnik,  Merrill  I.  Radar  Handbook.  New  York:  McGraw-Hill,  1970. 

38.  Skolnik,  Merrill  I.  Introduction  to  Radar  Systems.  New  York:  McGraw-Hill, 
1980. 

39.  Van  Trees,  Harry  L.  Detection,  Estimation,  and  Modulation  Theory,  Part  III. 
New  York:  John  Wiley  and  Sons,  1971. 

40.  Wilson,  David  J.  and  others.  “Laser  Atmospheric  Wind  Sounder  (LAWS)  Tech¬ 
nology  Status,”  Proceedings  of  SPIE  Fifth  Conference  on  Coherent  Laser  Radar: 
Technology  and  Applications,  1181:  79-83.  Bellingham,  WA:  SPIE,  1989. 

41.  Woodward,  Philip  M.  Probability  and  h formation  Theory,  with  Applications  to 
Radar.  New  York:  McGraw-Hill,  1955. 

42.  Wright  Research  and  Development  Center,  Air  Force  Systems  Command.  Con¬ 
tract  F33615-88-C-1756  with  Coherent  Technologies  Inc.  Wright-Patterson  AFB 
OH,  1  November  1989. 


96 


Vita 

Captain  Jerry  L.  Rogers  was  born  on  16  September  1959  in  Bowdon,  Georgia. 
He  graduated  from  Bowdon  High  School  in  1977  and  married  Rita  N.  Stewart.  He 
graduated  with  honors  from  West  Georgia  College  in  1982  earning  a  Bachelor  of 
Arts  degree  in  Mathematics.  After  teaching  high  school  mathematics  two  years, 
he  attended  the  USAF  officer  training  school  (OTS).  He  was  commissioned  in  May 
1984  and  served  his  first  tour  as  a  student  at  Auburn  University.  He  graduated  with 
honors  from  Auburn  in  March  1986  with  a  Bachelor  of  Science  degree  in  Electrical 
Engineering.  Upon  graduation,  he  served  three  years  at  Newark  AFB,  Ohio  as  a  test 
equipment  supervisor  and  inertial  navigation  systems  engineer.  He  was  selected  in 
March  1989  to  attend  the  School  of  Enginei.  ring.  Air  Force  Institute  of  Technology. 


Permanent  address:  2057  Bonnie  Da.le  Dr 
Bellbrook  OH  45305 


97 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  fcponiftg  burdcit  for  thit  «.ol^Ct1toni  gf  <nfcrirwn»cn  evumateg  to  i  fx*  «nt*udinq  t^e  tim»*  toi  ievn^win^j  injtruasoni.  >Cdi o  cn'Stm^  data  wurics. 

gathenn^and  wsmtammg  the  data  nfodcd,  and  4C<rpieting  ^jno  me  oi  int«,rmat<wn  >i»nd  v^.mTirnwtMatdmg  triisourden  oumate  v»  ♦ny  othui  a>pcs,t  vf  tbik 

COllfrCtiCn  of  »n|ormat»on,  irctudmg  tugQe>ticn»  to?  reduvtng  tnjy  4?grd**n  Wajf>»nqtun  Mt'aoqwar'Cf^  uiravioiaiefvt  -nfvtfnaitun  Oc*ctat«or*> and  ft  ‘pvrU.  IS  ^effcrwn 

Oavd Highway, Suite  l^tOa^Afhngtcn,  vA  ^2202‘4i0^^andtot»^eOttKe.jtManacj|cmt»niandi^wd-.jei>ipcrwQU»  (KcdMaJooPrvjett^0/C4-0t$8)^  Aashmgwn,  uC  aOS03 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 


December  1990 


3.  REPORT  TYPE  AND  DATES  COVERED 

_ Master’s  Thesis  _ 


4.  TITLE  AND  SUBTITLE 


Pulse  Compression  for  Aerosol  Ranging  With  Coherent  Pulse-Doppler  Lidar  Sys¬ 
tems 


6.  AUTHOR(S) 

Jerry  L.  Rogers,  Capt 


7.  PERFORMING  ORGANIZATION  NAM£(S}  AND  ADDRESS(ES} 

Air  Force  Institute  of  Technology,  WPAFB  OH  45433-6583 


9.  SPONSORING /MONITORING  AGENCY  NAME(S}  AND  AODRESS(ES) 

Electro-Optics 
WRDaAARI-2 
WPAFB  OH  45433-6543 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GE/ENG/90D-52 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  words) 

This  thesis  provides  a  thorough  review  of  conventional  matchedlfllter  radar  theory  as  it  applies  to  aerosol  sensing  lidar. 
Basic  matched-filter  radar  theory  and  the  complex,  dense  aerosol  target  model  eventually.lead  to  a  general  derivation  of 
the  matched-filter  radar  receiver  response  to  a  dense  aerosol  target  environment.  The  matched-filter  response  is  obtained 
in  terms  of  a  two  dimensional  convolution  of  the  target  scattering  function  and  the  radar  ambiguity  function.  The  range 
and  radial  velocity  resolution  of  various  radar  signals  were  compared  using  the  matched-filter  radar  receiver  and  scattering, 
function  models.  Pulse  con^ression  radar  signals  were  compart  to  the  simple  radar  pulses  cunently  in  use  with  existing 
lidar  systems,  and  in  each  case,  the  pulse  conpression  radar  signal  did  not  provide  a  significant  improvement  in  combined 
range  and  radial  velocity  resoluUon.  Tke^^ 


4.  SUBJECT  TERMS 


Pulse  Con^ression,  Aerosol  Ranging,  Coherent  Pulse4)oppler  Lidar,  Lidar  Wind  Sensing 


Systems,  bedaa^Laser  Radar 


Ranging,  com 


se^op] 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT  t 


Unclassified 


NSN  7540-01 -280-5500 


Unclassified 


Unclassified 


standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std  239- *8 
298-102 


